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I.  BACKGROUND 

The  U.S.  Navy  is  actively  engaged  in  promoting  the  application  of  state-of-the-art  and 
commercial-off-the-shelf  technologies  to  realize  cost-effective  platforms.  In  particular, 
research  and  development  is  focussed  on  improving  producibility,  enhancing  operational 
flexibility,  maximizing  survivability,  decreasing  manning  requirements  and  decreasing 
overall  system  cost.  From  the  power  system  perspective,  these  goals  have  translated  into 
research  initiatives  seeking  to  reduce  the  size  and  weight  of  the  principle  distribution 
components  and  to  optimize  the  operation  and  survivability  of  key  system  elements. 
Presently,  these  efforts  have  been  concentrated  in  four  areas:  zonal  distribution,  DC  Zonal 
Electric  Distribution  (DC  ZEDS),  Power  Electronic  Building  Blocks  (PEBBs),  and  the 
Integrated  Power  System  (IPS)  also  referred  to  as  the  Integrated  Electric  Drive  (IED). 

Zonal  Electrical  Distribution 

The  zonal  distribution  program  has  identified  that  a  zonal  architecture  has  a  number  of 
advantages  over  the  current  radial  distribution  architecture.  The  radial  network  includes 
generators  supplying  switchboards  then  a  myriad  of  feeder  cables  strewn  throughout  the 
ship  to  provide  power  to  vital  and  non-vital  electrical  loads.  The  zonal  architecture  is 
based  on  implementing  a  port  and  starboard  bus  and  sectioning  the  ship  service  electrical 
loads  into  a  number  of  zones  delineated  by  watertight  bulkhead  compartments.  Vital  loads 
within  a  zone  are  connected  to  either  bus  via  an  auctioneering  process.  The  advantages  of 
zonal  include  the  elimination  of  a  significant  amount  of  feeder  cables,  main  busses  only 
transition  watertight  compartments,  the  ship  may  be  fabricated  and  tested  in  zones,  and 
ship  construction  is  markedly  simplified.  The  arguments  for  zonal  distribution  are 
persuasive  and  the  only  issue  that  remains  is  whether  the  distribution  busses  are  three- 
phase  ac  or  dc. 

DC  Zonal  Electrical  Distribution 

The  DC  ZEDS  program  is  investigating  the  feasibility  of  zonal  dc  distribution.  In  DC 
ZEDS,  the  ac  generator  voltages  are  immediately  rectified  and  the  1200  Vdc  to  2000  Vdc 
is  sent  to  the  port  and  starboard  busses.  Each  bus  is  connected  to  an  electrical  zone 
through  a  power  converter  called  a  Ship  Service  Converter  Module  (SSCM).  The  SSCM 
buffers  the  main  bus  and  intra-zonal  loads,  monitors  zone  conditions,  and  adjusts  and 
regulates  the  main  bus  voltage  downward  to  a  level  commensurate  with  dc-to-ac  inverter 
requirements  (850-950  Vdc).  Three-phase  and  single-phase  ac  voltages  are  synthesized 
within  a  zone  by  a  power  converter  called  a  Ship  Service  Inverter  Module  (SSIM).  The 
SSIM  employs  "intelligent"  feedback  control  to  provide  tightly-regulated  ac  voltage  and 
current  to  the  corresponding  loads.  In  addition,  the  SSIM  on-board  intelligence  facilitates 
rapid  current  limiting  which  allows  for  a  degree  of  self  protection.  It  is  envisioned  that 
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standardized  converter  modules  will  be  paralleled  to  achieve  the  required  zone  power 
requirements  and  will  employ  soft-switching  technology  in  order  to  minimize  switching 
losses  and  maximize  control  bandwidth.  The  principle  benefits  of  DC  ZEDS  are  that  fault 
detection  and  clearing  are  both  simpler  and  faster  and  that  faults  are  now  isolated  to  a 
particular  zone.  The  time  lag  associated  with  initiating  ac  ' bus  transfers  is  virtually 
eliminated,  enhancing  the  integrity  of  power  flow  to  critical  electrical  apparatus.  In 
addition,  DC  ZEDS  eliminates  the  need  for  most  of  the  distribution  transformers  and  ac 
switchgear  providing  substantial  savings  in  size  and  weight.  The  ready  availability  of 
feedback-controlled  dc-to-ac  inverters  facilitates  the  application  of  variable  speed  drives 
which  in  turn  optimize  the  operation  of  blowers  and  pumps  while  implementing  an 
automatic  limiting  of  the  in-rush  currents  experienced  when  starting  large  motors.  DC 
ZEDS  also  eliminates  the  requirement  that  the  ship  service  generators  be  operated  at  60 
Hz.  This  allows  for  a  more  optimized  generator  design  in  terms  of  size,  weight  and  cost. 
Finally  SSCMs  and  SSIMs  are  multifunctional  and  allow  for  the  inclusion  of  more 
intelligent  power  management  and  fault  protection.  Thus,  despite  the  fact  that 
technological  hurdles  still  persist  in  realizing  a  dc  zonal  distribution  system,  the  advantages 
are  clear  and  the  arguments  are  sound. 

Power  Electronic  Building  Block 

The  PEBB  program  has  concentrated  on  developing  power  electronic  modules  which 
combine  a  power  section,  current  and/or  voltage  sensors  and  an  on-board  computer 
intelligence.  In  theory,  the  aggregate  can  be  programmed  and  interconnected  to  realize 
intelligent  circuit  breakers  or  more  complex  power  converters.  The  emphasis  thus  far  has 
been  on  using  existing  technology  to  realize  high-power  high-bandwidth  converters  for 
applications  such  as  DC  ZEDS.  The  high-bandwidth  control  is  realized  by  employing 
switching  frequencies  greater  than  20kHz  together  with  DSP-based  digital  controllers. 
Owing  to  the  intensive  input/output  (I/O)  requirements  of  many  of  the  DC  ZEDS 
components,  customized  DSP  resources  are  required  which  in  turn  require  algorithms  to  be 
developed  and  corresponding  code  to  be  generated.  The  development  of  reliable 
multifunctional  units  has  progressed  as  evidenced  by  several  testbed  demonstrations 
conducted  at  the  Naval  Surface  Warfare  Center  (NSWC)  in  Annapolis  Maryland. 

Integrated  Power  System 

Finally  the  IPS  program  has  revisited  the  feasibility  of  supplying  both  ship  service 
and  ship  propulsion  requirements  from  a  common  set  of  prime  movers.  A  current  surface 
combatant  utilizes  four  large  General  Electric  LM2500  gas  turbines  and  a  clutch  and 
reduction  gear  to  deliver  power  to  two  Controllable-Reversible-Pitch  (CRP)  propellers. 
Low  speeds  are  achieved  by  idling  the  turbine  and  adjusting  the  propeller  pitch.  Speeds 
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above  15  knots  are  achieved  by  a  combination  of  pitch  control  and  governed  throttle 
control  where  the  speed  of  the  turbine  is  adjusted  between  1200rpm  and  3600rpm.  Ship 
service  is  supplied  from  four  Allison  501K  gas  turbines  driving  2.5MW,  60Hz  synchronous 
machines.  Auxiliary  Propulsion  Units  (APUs)  are  available  to  provide  operation  up  to  5 
knots  and  are  powered  from  the  ship  service. 

In  electric  drive,  four  prime  movers  deliver  power  to  synchronous  machines  which  in 
turn  supply  both  ship  service  and  ship  propulsion  power  converters.  Propulsion  power 
converters  feed  propulsion  motors  which  directly  control  propeller  speed.  The  advantages 
of  IPS  are  manifold.  First,  the  number  of  prime  movers  is  reduced  which  leads  to 
associated  benefits  in  manning,  reliability,  cost  and  prime  mover  efficiency.  By  no  longer 
requiring  the  prime  mover  to  be  in  line  with  the  propeller  shaft,  the  gas  turbines  can  be 
located  on  higher  decks  reducing  the  size  of  inlet/outlet  ducting  and  minimizing  the 
probability  of  losing  power  during  a  flooding  casualty.  In  addition,  a  shorter  shaft  line 
reduces  the  complexity  of  aligning  equipment  across  several  ship  construction  zones. 
Since  flexible  electrical  connections  link  the  various  propulsion  apparatus,  any  prime  mover 
can  be  configured  to  power  any  propulsion  converter  and  either  propeller  shaft.  Also,  by 
eliminating  the  direct  connection  of  the  reduction  gears  and  gas  turbines  to  the  propeller 
shaft,  electric  drive  will  be  considerably  quieter  than  the  current  system,  reducing  ship 
signature.  From  a  fuel  consumption  point  of  view,  the  prime  mover  speed  is  decoupled 
from  the  propeller  shaft  speed  and  thus  the  prime  mover  may  be  operated  at  its  most 
efficient  setting  (around  3300rpm  for  an  LM2500).  This  in  turn  reduces  emissions.  The 
infinitely  variable  speed  control  available  with  electric  drive  makes  low-speed  maneuvering 
and  reversals  simple  and  rapid.  Finally,  IPS  also  provides  excess  generating  capacity  that 
may  be  used  for  powering  future  electrical  auxiliaries  such  as  Pulse-Energy  Weapons 
(PEWs)  and  the  ElectroMagnetic  Aircraft  Launching  System  (EMALS).  As  was  stated  for 
DC  ZEDS,  several  technical  hurdles  regarding  motor  and  converter  design  remain  to  be 
resolved;  however,  as  these  issues  are  confronted  by  R&D  resources,  BPS  will  be  available 
to  contribute  to  the  design  of  more  cost-effective  Naval  platforms. 

Naval  Postgraduate  School  Involvement 

During  the  past  six  years,  personnel  at  the  Naval  Postgraduate  School  (NPS)  have 
supported  the  research  efforts  into  zonal  distribution,  DC  ZEDS,  PEBB  and  BPS  through  a 
number  of  reimbursable  research  projects  and  over  twenty  Master’s  theses.  Currently,  the 
principle  focus  of  the  faculty  has  been  on  generating  deliverables  concerning  DC  ZEDS 
and  PEBB.  In  particular,  NPS  faculty  have  supervised  the  fabrication  of  several  medium- 
power  SSCM  power  sections,  the  design  and  validation  of  a  number  of  SSCM  feedback 
algorithms,  the  implementation  of  SSCM  regulators  using  both  analog  and  DSP  resources, 
the  analysis  and  implementation  of  various  paralleling  regimes  for  SSCMs,  the 
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development  of  closed-loop  algorithms  for  SSIMs,  the  study  of  the  interconnection 
feasibility  and  dynamics  of  multiple  SSCMs  and  SSIMs,  the  analysis  and  fabrication  of  a 
soft-switching  SSCM,  the  analysis  and  simulation  of  a  soft-switching  SSIM,  and  the 
investigation  of  the  effects  of  input  filter  parameters  on  the  stability  of  the  dc  distribution 
system.  A  number  of  conference  and  journal  publications  have  resulted  from  these 
research  efforts. 

In  terms  of  student  participation,  the  above  efforts  have  been  extensively  supported  by 
over  twenty  thesis  projects.  Kline  (Dec.  1993)  developed  a  detailed  SIMULINK  model  of 
a  representative  portion  of  a  shipboard  electric  power  distribution  system.  He  addresesed 
the  modeling  and  control  of  a  three-phase  synchronous  generator.  In  Colby  (Dec.  1993), 
the  constant  power  characteristics  of  DC  ZEDS  was  investigated  and  some  observations 
regarding  stability  and  controllability  were  made  upon  employing  a  reduced-order 
representation  of  the  dc-dc  converter  modules  using  PSPICE.  Mark  Kipps  (Mar.  1994) 
developed  a  means  of  using  a  built-in  differential  algebraic  solver  within  ACSL  to  simulate 
stiffly-connected  isolated  power  systems  without  relying  on  reformulated  machine 
representations  or  introducing  fictitious  circuit  components.  Blalock  (Mar.  1995)  initiated 
some  of  the  detailed  modeling  of  a  DC  ZEDS  system  consisting  of  a  phase-controlled 
rectifier,  dc-dc  buck  chopper  and  three-phase  inverter.  In  addition,  he  performed  the 
preliminary  investigations  into  various  closed-loop  algorithms  for  the  buck  chopper  and 
implemented  a  bench-top  hardware  testbed.  He  then  incorporated  feedback  control  of  the 
dc-dc  converter  using  the  dSPACE  hardware-in-the-loop  development  system.  Blalock 
further  conducted  hardware  and  software  studies  and  used  the  results  to  validate  his 
computer  representations. 

Paralleling  and  interleaving  dc-dc  buck  choppers  was  first  considered  by  Filor  (Sep. 
1995)  in  a  detailed  ACSL  simulation  exercise  that  included  a  model  of  a  steam  turbine 
driven  synchronous  machine,  rectifier,  filter  and  buck  choppers.  Some  of  the  issues 
regarding  paralleing  were  brought  to  the  fore  for  treatment  by  future  students.  Salerno 
(Jun.  1996)  significantly  advanced  the  proposed  closed-loop  dc-dc  converter  algorithm, 
integrated  a  DSP  solution,  coded  the  algorithms,  aided  in  implementing  the  hardware  setup 
and  conducted  a  series  of  validation  studies.  Close  to  the  same  time,  Oberley  (Dec.  1996) 
developed  a  detailed  ACSL  representation  of  an  ARCP  inverter.  Then  using  reduced-order 
representations,  investigated  various  stationary  and  synchronous  reference  frame  current 
control  algorithms  using  a  hardware  testbed  and  the  dSPACE  hardware-in-the-loop 
development  system.'  In  addition,  he  investigated  and  contrasted  the  PWM 
implementations  using  Sine-Triangle  PWM  (STPWM)  and  Space  Vector  Modulation 
(SVM).  These  studies  were  conducted  in  ACSL  assuming  an  RL  load.  Finally,  Oberley 
considered  an  involved  feedback/feedforward  voltage  control  scheme  which  involved 
significant  reference  frame  theory.  In  Nelson  (Dec.  1996),  the  one-cycle  control  algorithm 
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for  a  buck  chopper  was  considered  and  implemented.  Comparisons  were  made  between 
the  computer  representation  and  the  hardware.  Allen  (Mar.  1997)  documented  the  design 
and  fabrication  of  several  dc-dc  buck  chopper  power  sections  while  Badorf  (Jun.  1997) 
detailed  the  controller  design,  implemented  the  required  gating  circuitry,  and  built  and 
documented  the  associated  analog  control  hardware.  Badorf  also  developed  ACSL 
simulations  of  a  variety  of  interconnection  possibilities  for  DC  ZEDS  and  performed 
various  transient  studies. 

Langlois  (Jun.  1997)  supervised  the  fabrication  of  a  PEBB  testbed  used  for 
interconnecting  buck  choppers  and  ARCP  inverters.  He  then  compiled  some  hardware 
studies  investigating  the  transient  response  of  the  units  for  different  configurations  and 
loadings.  In  Hanson  (Jun.  1997),  the  operation  of  the  PUC  was  documented  and  the 
closed-loop  algorithm  for  the  buck  chopper  was  programmed  and  validated.  In  addition, 
Hanson  set  forth  the  approach  for  implementing  the  closed-loop  control  of  the  ARCP 
inverters.  Floodeen  (Sep.  1998)  further  delved  into  the  operation  of  the  PUC  and 
investigated  developing  C++  programs  for  controlling  the  various  converters.  In  addition, 
Floodeen  incorporated  several  auxiliary  functions  into  the  buck  chopper  control,  including 
an  non-droop  based  algorithm.  He  also  set  forth  the  method  and  coding  required  to 
implement  the  closed-loop  current  control  of  the  ARCP.  Presently,  Moore  (Mar.  1999)  is 
investigating  a  de-centralized  buck  chopper  paralleling  algorithm  that  relies  on  injected 
signal  detection.  He  plans  to  develop  a  hardware  testbed  to  illustrate  proof-of-concept. 
Marinac  (Sep.  1999)  will  extend  Moore’s  results  to  the  paralleling  of  ARCP  inverters. 
Greseth  (Sep.  1999)  is  investigating  the  stability  issues  regarding  the  dc-dc  converter  input 
filters  and  generator/rectifier  source.  In  particular,  the  stability  problems  associated  with 
the  constant  power  loads  is  analyzed  and  simulated  using  ACSL.  A  hardware  testbed  is 
used  to  validate  conclusions. 

There  have  been  a  host  of  other  thesis  topics  including  induction  motor  vector  control 
(Pierce,  Sep.  1995),  cycloconverter  drives  for  ship  propulsion  (Mercer,  Dec.  1996)  and  a 
paper  study  on  integrated  electric  drive  design  issues  (Arrington,  Sep.  1998)  that  have 
contributed  to  various  aspects  of  IPS.  Work  has  been  performed  on  the  simulation  and 
control  of  solar-powered  vehicles  (Roerig,  Mar.  1995  and  Yourkowski,  Mar.  1996)  and  the 
design  and  development  of  a  UPS  for  a  marine  vehicle  (Garcia,  Dec.  1994  and  Callahan, 
Mar.  1997).  Zupfer  (Dec.  1993)  and  Jones  (Jun.  1995)  investigated  Active  Power  Line 
conditioner  issues  while  Hand  (Mar.  1994)  was  involved  in  the  development  of  a  power 
supply  for  an  amateur  satellite.  Finally,  hardware  setups  were  used  to  perform 
investigations  into  a  slip  energy  recovery  system  (Tait,  Dec.  1995)  and  a  Mapham 
converter  (Tyner,  Dec.  1995). 
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The  remainder  of  this  report  focusses  on  the  deliverables  pertaining  to  a  control  law 
library  for  PEBB  devices.  The  following  sections  contain  detailed  descriptions  of  closed- 
loop  algorithms  for  dc-dc  converters,  dc-ac  inverters  and  ac-dc  phase-controlled  rectifiers. 
These  algorithm  are  designed  to  be  realizable  with  the  PEBB  Universal  Controller 
(Generation  0)  and  have  been  validated  using  detailed  and  reduced-order  Advanced 
Continuous  Simulation  Language  (ACSL)  simulation  models. 
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A.  ARCP  Inverter  Control 

The  Auxiliary  Resonant  Commutated  Pole  (ARCP)  inverter  operates  similar  to  a 
classical  three-phase  bridge  inverter  except  that  at  the  switching  instances,  auxiliary 
circuitry  ensures  that  each  semiconductor  device  is  activated  with  approximately  zero  volts 
across  it.  For  a  given  output  power  level,  the  transistor  soft-switching  facilitates  higher 
switching  frequencies.  By  extending  switching  frequencies  from  about  20kHz  up  to  60kHz 
for  units  rated  in  the  tens  of  kilowatts,  the  ARCP  inverter  allows  for  faster-acting  higher- 
bandwidth  feedback  control  systems  to  regulate  the  voltages  across  or  currents  through  a 
three-phase  load  and  for  smaller  output  filter  elements. 

An  ARCP  inverter  consists  of  six  main  transistors  with  freewheeling  diodes  and  a  set 
of  six  auxiliary  transistors  with  corresponding  freewheeling  diodes.  For  the  set  of  units 
delivered  to  NSWC  (via  NPS)  by  Penn  State  University  (PSU),  the  user  has  the  option  of 
either  generating  all  twelve  device  gate  signals  or  specifying  only  the  main  transistor  gate 
signals  and  allowing  the  ARCP  unit  to  appropriately  activate  the  auxiliary  devices.  All 
actuation  signals  are  optically  coupled  to  the  ARCP  to  maximize  isolation  and  optimize 
noise  immunity.  The  Power  Electronic  Building  Block  (PEBB)  Universal  Controller 
Generation  0  (PUC-GENO)  was  developed  with  extensive  Input/Output  (I/O)  capability  (12 
optical  output  transmitters  and  10  Analog-to-Digital  converter  input  ports)  to  interface  with 
the  ARCP  units.  Currently,  the  PUC-GENO  is  Assembly  Language  programmable'  with 
future  generations  anticipated  to  be  programmable  in  C.  In  the  remainder  of  this  section,  it 
is  assumed  that  the  auxiliary  device  control  was  left  to  the  on-board  ARCP  circuitry  so  that 
the  PUC-GENO  can  be  dedicated  to  implement  more  sophisticated  closed-loop  algorithms. 
This  seems  reasonable  since  future  inverter  algorithms  may  also  require  paralleling  and 
condition  monitoring  features  that  would  exhaust  the  I/O  capability  of  the  PUC-GENO  if  it 
was  also  required  to  control  the  auxiliary  switches. 

Closed-loop  inverter  control  is  a  mature  subject  area.  The  intent  of  this  section  is  to 
explore  a  number  of  closed-loop  voltage  and  current  algorithms  that  would  be  appropriate 
for  three-phase  naval  shipboard  loads.  First,  we  will  focus  on  how  a  duty  cycle  signal  may 
be  derived  using  Sine-Triangle  Pulse-Width-Modulation  (STPWM)  and  then  using  Space 
Vector  Modulation  (SVM).  A  brief  exposition  on  reference  frame  theory  will  be  included 
with  SVM.  Reference  frame  analysis  is  also  key  to  the  current  and  voltage  control 
stratagems  presented  in  the  subsequent  sections.  The  current  control  algorithms  are 
addressed  first,  documenting  the  control  in  the  stationary  reference  frame  then  extending 
the  analysis  to  the  synchronous  reference  frame.  The  voltage  control  algorithm  is  a 
sophisticated  control  loop  incorporating  both  feedback  and  feedforward  elements.  Both  the 
voltage  and  current  controls  can  be  ultimately  implemented  using  either  STPWM  or  SVM. 
Finally,  two  algorithms  for  implementing  a  current-limited  induction  machine  startup  are 
documented. 
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1.  Modulation  Strategies 

The  most  basic  inverter  modulation  strategies  are  the  180°  (six-step)  and  120°  modes. 
The  180°  mode  involves  gating  the  inverter  leg  switches  in  a  complementary  fashion  with 
each  device  on  for  180°  of  an  output  cycle.  The  signals  for  the  three  inverter  legs  are 
phase  displaced  by  120°  to  achieve  a  balanced  set.  The  120°  mode  gates  each  device  in  an 
inverter  leg  for  120°  of  an  output  cycle  with  60°  separating  the  gating  of  the  upper  and 
lower  switches.  Once  again,  the  signals  for  the  three  inverter  legs  are  phase  displaced  by 
120°  to  achieve  a  balanced  set.  Both  of  these  algorithms  are  simple  but  are  unattractive 
since  the  fundamental  output  amplitude  is  not  directly  controllable  and  considerable 
harmonics  are  injected  into  the  load.  The  solution  is  to  intelligently  place  notches  into 
these  quasi-rectangular  waveforms  using  one  of  the  variants  of  Pulse-Width-Modulation 
(PWM). 


a.  Sine-Triangle  Pulse-Width-Modulation 

STPWM  employs  a  high-frequency  carrier  waveform  together  with  a  balanced  set  of 
modulating  signals  at  the  desired  output  frequency.  Each  modulating  signal  controls  the 
transistor  switching  for  a  given  leg.  When  the  modulating  signal  is  greater  (less)  than  the 
carrier  waveform,  the  upper  (lower)  switch  is  gated.  Three  modes  of  operation  are  possible 
depending  on  the  relative  amplitudes  of  the  carrier  and  modulation  signals.  If  the 
modulating  signal  amplitudes  are  less  than  the  earner  amplitude,  STPWM  is  operating  in 
the  linear  mode  and  the  phase  voltage  amplitude  is  directly  proportional  to  the  modulating 
signal  amplitude.  For  instance  if  Vpk  denotes  the  peak  value  of  the  modulating  signal,  At 
the  peak  triangle  amplitude  and  Vdc  the  inverter  dc  input  voltage,  the  fundamental  phase 
voltage  is  given  by 

Vas’!  =  VpkI^cos((0t)  .  ^ 

where  a>  is  the  desired  output  frequency  in  radians/sec.  When  the  modulating  signal 
amplitude  exceeds  the  triangle  amplitude,  STPWM  transitions  into  the  pulse-dropping 
(over-modulation)  mode  where  there  is  a  nonlinear  relationship  between  the  modulation 
amplitude  and  the  phase  voltage  amplitude.  If  we  define  the  modulation  index 


M  = 


Vpk 

A* 


(2) 


then 


vas,l 


VdcM 


Sin_,(l^)  +  iT'  (M)2 


cos(  CO  t ) 


(3) 


M  7  '  M  y '  '  M 

Finally,  if  the  modulating  amplitude  is  significantly  larger  than  the  carrier  amplitude,  only 
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two  waveform  intersections  occur  per  desired  output  cycle  and  STPWM  operates  in  180° 
mode  and 


vas,l 


2Vdc 

n 


COS(O)t) 


(4) 


In  any  of  the  modes  as  long  as  the  modulating  signals  are  not  corrupted  by  harmonics  and 
are  at  a  much  lower  frequency  than  the  carrier  waveform,  a  given  inverter  transistor  is 
switched  "on"  and  "off  once  per  cycle  of  the  carrier  waveform.  This  fixes  the  switching 
frequency  as  compared  to  a  hysteresis  inverter  control  and  also  produces  well-defined 
harmonics. 

Normally,  the  modulating  signals  used  in  STPWM  are  derived  from  feedback  and 
possibly  feedforward  signals.  The  strategy,  together  with  protections  against  gating  both 
transistors  in  a  given  inverter  leg,  is  readily  implemented  using  standard  Integrated  Circuit 
(IC)  chips.  However,  it  is  intended  that  the  PUC-GENO  be  used  to  control  the  ARCP  so  it 
is  necessary  to  discretize  the  STPWM  algorithm. 

If  we  assume  that  we  calculate  the  modulating  signals  each  time  we  sample  the 
feedback  quantities,  we  implicitly  assume  that  that  modulating  signal  stays  constant  over  a 
switching  (triangle)  cycle  and  the  inverter  leg  duty  cycle  can  be  evaluated  as  follows. 
Designating  D^,  DtB,  and  DtC  as  the  duty  cycle  amounts  for  each  of  the  top  switches  in  the 
inverter  legs, 


DtA 

DtB 

Dtc 


A, 

"t”  Vmod,a 

2  A, 

A, 

Vmod,b 

2  A, 

A, 

Vmod,c 

2  A, 

(5) 

(6) 
(7) 


Theoretically  the  duty  cycles  are  constrained  to  lie  between  zero  and  one;  however,  the 
turn-off  and  turn-on  times  of  the  transistors  disqualify  duty  cycles  that  result  in  very 
narrow  pulses.  If  is  the  minimum  time  length  achievable  for  a  pulse,  then 


Dmin 


tmin 

T 

2 


(8) 


where  y  is  half  of  a  switching  period  and 

Dmax  =  1  -  D™  (9) 

Thus,  if  a  duty  cycle  falls  betweem  and  one,  we  need  to  round  the  desired  value  up  to 
one,  effectively  eliminating  that  transition  of  the  modulating  and  carrier  waveforms.  If  a 
duty  cycle  falls  betweem  D^n  and  zero,  we  need  to  round  the  desired  value  down  to  zero. 
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The  duty  cycles  can  then  be  converted  into  counts  for  the  PUC-GENO  counter/timers 
and  implemented  as  described  by  Floodeen.  Issues  that  need  to  be  considered  are  (1)  the 
quantization  effect  of  converting  the  duty  cycle  to  a  count,  (2)  the  harmonic  performance 
of  the  inverter  in  pulse-dropping  mode,  (3)  the  effect  of  the  delay  inserted  between  gating 
"off"  a  leg  switch  and  gating  "on"  the  other  leg  switch,  (4)  the  selection  of  At  together  with 
the  feedback  gains,  and  (5)  the  effect  on  harmonics  of  using  asynchronous  STPWM  where 
the  switching  frequency  is  not  a  multiple  of  the  desired  inverter  output  frequency. 

b.  Space  Vector  Modulation 

A  three-phase  bridge  inverter,  as  shown  in  Figure  1,  has  eight  possible  switching 
states  as  illustrated  in  Table  1  where  it  is  assumed  that  the  leg  transistors  are  fired  in  a 
complementary  fashion.  Also,  a  ’1’  in  the  table  denotes  that  the  switch  is  closed;  a  ’O’  in 
the  table  denotes  that  the  switch  is  open. 

Table  1.  Possible  Inverter  Switching  States 


State 

Si 

s2 

S3 

1 

1 

0 

0 

2 

1 

1 

0 

3 

0 

1 

0 

4 

0 

1 

1 

5 

0 

0 

1 

6 

1 

0 

1 

7 

1 

1 

1 

8 

0 

0 

0 

The  voltages  across  the  lower  transistor  switches  (vap,  vbp  and  v^,)  are  known  once  the 
gating  signals  are  specified.  That  is,  if  the  lower  switch  is  gated  vap  =  0,  while  if  the  upper 
switch  is  gated  vap  =  Vdc.  It  can  be  shown  that  for  balanced  loads,  including  symmetrical 
induction  machines  and  synchronous  machines,  it  is  always  true  that 

Vas  +  Vbs  +  Vcs  =  0  (10) 

which  then  allows  us  to  formulate  the  phase  voltages  based  solely  on  knowledge  of  vap,  vbp 


and  vcp. 


vas  _  g  VaP  3  Vbp 


(ID 
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(12) 


vbs  -  j  VfcP  3  VaP  3  v<=p 

VCS  =  J  Vcp  -  J  vap  —  —  Vbp  (13) 

The  phase  voltages  for  each  inverter  state  can  be  found  by  substitution  as  illustrated  in 
Table  2. 


Table  2.  Phase  Voltages  for  Each  Inverter  State 


State 

vap 

Vbp 

Vcp 

Vas 

Vbs 

Vcs 

1 

Vdc 

0 

0 

2Vdc 

3 

Vdc 

3 

Vdc 

3 

2 

Vdc 

Vdc 

0 

Vdc 

3 

Vdc 

3 

2Vdc 

3 

3 

0 

Vdc 

0 

Vdc 

3 

2  Vdc 

3 

Vdc 

3 

4 

0 

Vdc 

Vdc 

2Vdc 

3 

Vdc, 

3 

Vdc 

3 

5 

0 

0 

Vdc 

Vcfc 

3 

Vdc 

3 

2Vdc 

3 

6 

Vdc 

0 

Vdc 

3 

2Vdc 

3 

Vdc 

3 

7 

Vdc 

Vdc 

Vdc 

0 

0 

0 

8 

0 

0 

0 

0 

0 

0 

The  six-step  algorithm  is  implemented  by  sequentially  selecting  the  first  six  states,  each 
assigned  for  \  of  the  desired  output  period. 

Reference  frame  theory  is  frequently  applied  in  the  analysis  of  electric  machinery  in 
order  to  achieve  marked  simplification  of  the  modeling  equations.  The  transformations  or 
change  of  variables  are  diffeomorphic,  one-to-one,  and  simply  allows  us  to  view  the 
machine  dynamics  in  a  different  coordinate  system.  The  stationary  reference  frame  is 
frequently  used  in  the  analysis  of  induction  machines.  It  consists  of  an  algebraic 
relationship  where  fictitious  q  and  d-quantities  are  derived  from  the  actual  phase  quantities. 
The  phase  voltage  quantities  are  transformed  into  the  stationary  reference  frame  by 


Vqs  =  Va 


V3 


Vbs  + 


V3 


(14) 


(15) 
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Figure  1.  Three-Phase  Voltage  Source  Inverter 


If  we  view  the  stationary  q-axis  as  being  a  horizontal  axis  and  the  stationary  d-axis  as 
being  orthogonal  and  directed  downward,  we  may  then  characterize  the  instantaneous 
values  of  Vq*  and  vjs  by  a  vector  in  the  qd-plane  where  the  amplitude  of  that  vector  is  given 

by 


(16) 


while  the  angle  is  given  by 


^qAang  —  tan  (  «  ) 


Therefore,  each  inverter  switching  state  corresponds  to  unique  values  of  vq*  and  v*  and  can 
thus  be  represented  by  a  unique  vector  in  the  qd-plane.  These  are  termed  "space  vectors" 
and  are  referred  to  herein  as  V!  through  V,.  The  vector  descriptions  are  listed  in  Table  3 
and  are  illustrated  in  Figure  2.  As  illustrated,  the  first  six  inverter  space  vectors  are  equal 
amplitude  and  are  evenly  displaced  in  the  qd-plane.  States  7  and  8  are  represented  by  zero 
amplitude  vectors. 

The  SVM  idea  stems  from  exploiting  the  inverter  switching  state  representation  in  the 
qd-plane.  If  we  can  generate  desired  values  for  vq*  and  v|,  then  we  can  implement  those 
values  by  appropriately  averaging  the  various  inverter  space  vectors.  Let’s  define  Tsw>per  as 
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Table  3.  Space  Vector  Representations  of  Inverter  States 


State 

v s 

vqs 

Vds 

^qd,amp 

^qd,ang 

1 

2Vdc 

3 

0 

2  Vdc 

3 

0° 

2 

Vdc 

3 

Vdc 

€ 

2  Vdc 

3 

60° 

3 

Vdc 

3 

Vdc 

V3 

2Vdc 

3 

120° 

4 

2Vdc 

3 

0 

2  Vdc 

3 

180° 

5 

vdc 

3 

Vdc 

V3 

2  Vdc 

3 

240° 

6 

Vdc 

3 

Vdc 

V3 

2  Vdc 

3 

300° 

7 

0 

0 

0 

origin 

8 

0 

0 

0 

origin 

the  switching  period.  It  will  correspond  to  the  amount  of  time  allocated  to  the  various 
inverter  states  used  to  synthesize  the  desired  space  vector.  The  qd-plane  is  divided  into  six 
sectors  defined  by  v^ang-  In  order  to  minimize  Tswper  relative  to  the  desired  output  period, 
we  choose  to  only  use  states  adjacent  to  the  sector  in  which  the  desired  state  vector 
resides.  Also,  a  sequence  of  states  exists  which  minimizes  the  number  of  switchings  per 
cycle.  These  are  illustrated  in  Table  4.  Note  that  the  optimal  switch  sequence  uses  both 
zero  vectors,  States  7  and  8.  The  averaging  process  may  be  viewed  as  follows.  The 
weighting  of  adjacent  inverter  states  is  used  to  achieve  the  phase  (y)  of  the  desired  space 
vector.  The  weighting  of  these  states  relative  to  the  zero  states  is  used  to  achieve  the 
desired  space  vector  amplitude.  Let’s  consider  Sector  I  and  define 

Tsw,per  =  T]  +  T2  +  Tz  (18) 

where  T,  is  the  time  spent  at  State  1,  T2  is  the  time  spent  at  State  2,  and  Tz  is  the  time 
spent  at  the  zero  state.  Note,  the  sequences  illustrated  in  Table  4  consist  of  two  such 

switching  periods  --  one  starting  at  State  8  and  one  starting  at  State  7.  If  we  next  define  V, 

as  the  space  vector  for  State  1  and  V2  as  the  space  vector  for  State  2,  then  the  average 
space  vector  (V)  achieved  over  a  switching  period  is  given  by 

T,Vi  +  T2V2 

T 

A  sw,per 


V 


(19) 
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State  3  State  2 


Stationary 

q-axis 


Figure  2.  Illustration  of  Space  Vectors 


Table  4.  Optimum  Switch  Sequences 


Sector 

^qd,ang 

Switch  Sequence 

I 

0°  <  y  <  60° 

8, 1,2, 7,2,1,... 

n 

60°  <  y  <  120° 

8, 3, 2, 7, 2, 3,... 

m 

120°  <  y  <  180° 

8, 3, 4,7,4, 3,... 

IV 

180°  <  y  <  240° 

8, 5, 4, 7,4, 5,... 

V 

240°  <  y  <  300° 

8, 5, 6, 7,6, 5,... 

VI 

300°  <  y  <  360° 

8, 1,6, 7,6,1,... 

The  assembly  of  the  average  space  vector  is  illustrated  in  Figure  3.  This  expression  is  set 
equal  to  the  desired  space  vector  V*  which  can  be  decomposed  into  qd-components  (in 
Sector  I)  as 

vqs*  =  Vpkcosy  (20) 

y£  =  -VpksinY  •  (21) 

where  y  is  the  angle  of  the  desired  space  vector  (equation  17).  Substituting  into  the  qd- 
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descriptions  for  V,  and  V2  (Table  3)  yields  two  equations  in  the  two  unknowns  Tj  and  T2. 


T,  =  T. 


V, 


sw.per  y 


3  V3  . 

-cosy  - 


•  vpk  r- 

T2  =  TSW(Per  -f-  V3  siny 


f  dc 


and  then  the  zero  state  time  is  found  from 
Tz  =  Tswper  -  Tj  -  T2 


(22) 

(23) 

(24) 


Figure  3.  Derivation  of  the  Average  Space  Vector 


Given  that  we  define 


t  =  T 

La  Asw,per  y 


qs 


dc 


Vds* 


Tp  -  Tsw  per  y 


dc 


(25) 

(26) 


we  can  establish  the  switch  times  for  the  devices  in  each  sector.  These  are  summarized  in 
Table  5.  As  mentioned  above,  the  optimal  sequence  for  a  given  sector  uses  both  zero  states 
as  illustrated  in  Table  4.  For  reasons  of  symmetry,  half  of  the  zero-state  time  illustrated  in 
Table  5  is  allocated  to  State  8  and  half  to  State  7.  Thus,  the  appropriate  sequence  may 
begin  with  either  zero  state.  The  corresponding  duty  cycles  for  the  top  switches  in  each 
leg  are  documented  in  Table  6.  Whether  the  duty  cycle  is  initially  high  or  low  for  a  given 
cycle  depends  on  the  current  zero  state.  That  is,  if  the  previous  sequence  ends  on  State  8, 
then  the  upper  switch  remains  low  (open)  for  l  -  D  of  the  cycle.  If  the  sequence  ends  on 
State  7,  the  upper  switch  remains  high  (closed)  for  D  of  the  cycle. 

One  issue  that  remains  in  implementing  SVM  is  to  establish  what  happens  if  the 
desired  space  vector  has  an  amplitude  larger  than  that  achievable.  This  corresponds  to  the 
over-modulation  region  for  STPWM.  For  Sector  I,  a  given  Vpk  is  realizable  as  long  as  the 
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Table  5.  Switch  Time  Calculations 


Sector 

Non-Zero  State  Times 

Zero-State  Time 

I 

3  V3 

Tl  =  1 ta  +  Ttp 

T2  =  -V3tp 

TSw,per  -  T,  -  T2 

n 

T  3t  ^  t 

T2  -  -  2  tp 

3,  V3 

T3  -  2  ta  2  p 

Tsw,per-T2-T3 

in 

T3  = 

3  V3 

t4  =  -■ jia  +  -ytp 

TSw,per  T3  —  T4 

TV 

3  V3  t 

T4  -  -  2  2  p 

T5  =  V3tp 

Tsw,per  -  T4  -  T5 

Y 

3  V3  t 

Ts  =  -y  la  +  -ytp 

3  V5 

T6  =  +  -ytp 

Tsw,per-T5-T6 

VI 

t6  =  ^tp 

„  3  V3 

Tl  ~  2ta  "  2  15 

Tsw,per  ~  T6  -  T] 

following  condition  is  satisfied 


vpk  < 


Vdc 

V3 


1 

cos(y  -  30°) 


(27) 


implying  that  a  sinusoid  of  ~  can  be  perfectly  tracked  through  a  sector  (14%  larger  than 

in  STPWM).  If  (27)  is  violated,  inadmissible  switch  times  will  result.  Two  options  exist: 
the  desired  space  vector  can  be  scaled  down  to  the  maximum  amplitude  given  in  (27), 
retaining  the  desired  phase  relationship,  or  the  closest  inverter  state  can  be  engaged  for  the 
entire  Tswper.  The  determination  can  be  based  on  which  solution  results  in  the  smaller  error 
vector.  For  instance  using  the  scaling  approach  for  Sector  I,  we  would  have  the  maximum 
vector  amplitude  constrained  to  be 


VLim  V3  cos(y  -  30°) 

and  therefore  the  scaled  desired  qd-voltages  would  be  given  by 


(28) 
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Table  6.  Top  Switch  Duty  Cycles 


Sector 

Dtfi 

Dtc 

I 

T, 

Tz 

+  X2  +  T 

Tz 

Tj+t 

Tz 

2 

T 

A  sw,per 

T 

A  sw.per 

T 

A  sw,per 

n 

Tz 

t’  +  T 

Tz 

T2  +  T3  +  — 

T 

Az 

T 

T 

A  sw,per 

T 

A  sw,per 

T 

A  sw,per 

m 

T 

Az 

T 

Tz 

T3  +  T4  +  Y 

Tz 

t,  +  T 

T 

A  sw,per 

Tsw,per 

Tsw,per 

IV 

T. 

2 

Tz 

t4+t 

Tz 

T4  +  T5  +  -f 

T 

A  sw,per 

T 

A  sw,per 

T 

A  sw,per 

V 

Tz 

T6+T 

Tz 

2 

Tz 

Ts  +  T6  +  Y 

T 

A  sw,per 

T 

A  sw,per 

T 

A  sw,per 

VI 

Tz 

+  Te  +  T 

t2 

2 

Tz 

t‘  +  T 

T 

A  sw.per 

T 

A  swtper 

T 

A  sw,per 

Vqs,new  ~  Vym  COS  Y  (29) 

vds,new  =  VLim  Y  (3(^) 

these  would  then  be  used  to  calculate  switch  times  Ti  and  T2.  The  resulting  error 
amplitude  would  be  given  by 

Veir.ampl  =  ^pk  —  ^Lim  (31) 

Expressions  for  VUm  in  the  remaining  sectors  are  documented  in  Table  7.  Assuming  y  lies 
between  -30°  and  30°  where  State  1  is  the  closest,  the  voltage  amplitude  error  in  this  case 
would  be 

VOT,amp2  =  a/(^)2  +  (Vp*k)2  -  j  vdc  v;k  cosy  (32) 

Expressions  governing  this  voltage  error  for  the  remaining  values  of  y  are  listed  in  Table  8. 
Equations  (31)  and  (32)  can  be  appropriately  compared  to  assess  which  solution  minimizes 
the  voltage  amplitude  error.  The  SVM  algorithm  may  viewed  in  the  following  steps: 
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1.  The  control  specifies  commanded  values  of  vqs*  and  vqs* 

VS* 

2.  From  which  we  can  evaluate  y  =  tan'1  ( ),  the  angle  of  the 
desired  space  vector. 

3.  Determine  the  sector  based  on  y  (Table  4). 

4.  Determine  the  desired  peak  value  of  the  space  vector 

vPk  =  V(vqs’)2  +( vaT ) 2 

5.  Calculate  the  maximum  achievable  space  vector  for  the  y  (VLim) 
from  Table  7. 

6a.  If  VPk  is  less  than  then  calculate  ta  and  tp  as  directed  in  (25) 
and  (26). 

6b.  If  VP1  is  greater  than  VUm  then  calculate  scaled  versions  of  the 
qd-commanded  voltages  from  (29)  and  (30)  which  are  then  used  in 
(25)  and  (26)  to  establish  ta  and  tp.  (Unless  the  error  vector  using  a 
single  space  vector  "full  on"  is  smaller.  Then,  set  the  corresponding 
duty  cycles  to  one  or  zero  for  the  given  leg  and  skip  the  remaining 
steps.) 

7.  Identify  the  switch  times  from  Table  5. 

8.  Identify  the  switch  leg  duty  cycles  from  Table  6. 

9.  Round  duty  cycles  <  D  <  1  up  to  one  and  duty  cycles 
0  <  D  <  Dmin  down  to  zero. 

10.  If  presently  at  State  8,  the  upper  switch  stays  low  (off)  for  1  -  D. 
If  presently  at  State  7,  the  upper  switch  stays  high  (on)  for  D. 


11.  Repeat  each  TSWiPer. 
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Table  7.  Maximum  Achievable  Space  Vector  for  Each  Sector 


Sector 

Vum 

I 

Vdc 

1 

V3 

cos(y-  30°) 

n 

Vdc 

1 

V3 

cos(y-  90°) 

m 

Vdc 

i 

"V3 

cos(y-  150°) 

IV 

Vdc 

i 

V3 

cos(  y  -  210° ) 

V 

Vdc 

1 

V3 

cos(  7-270°) 

VI 

Vdc 

1 

V3 

cos(y-  330°) 

Table  8.  Error  Vector  Amplitudes  Using  Single  State 


Desired  Vector  Location 

Verr,amp2 

-30°  <  y  <  30° 

30°  <  Y  <  90° 

90°  <  Y  <  150° 

150°  <  Y  <  210° 

210°  <  Y  <  270° 

270°  <  Y  <  330° 

3  )2  +  V  3  Vpkcos(y) 

A/(^)2  +  Vp*k2  -  3dc  Vpkcos(y  60°) 

+  Vp*k2  -  — Vpkcos( Y -  120° ) 

V(^>2  +  Vpk2  -  ^v;kcos(Y-180°) 

+-V£  -  — " Vpt cos( Y - 240° ) 

( -f-  f  +  Vp*k2  -  3  dc  Vp;  cos(  Y  300° ) 

c.  Hysteresis  ( Bang-Bang )  Modulation 

Hysteresis  or  Bang-Bang  PWM  is  implemented  by  having  the  error  between  the 
commanded  abc-quantities  and  the  actual  abc-quantities  directly  dictate  the  switch  status  in 
each  inverter  leg.  For  instance  for  Leg  A, 
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If  C  -  ias  >  £>  then  close  the  top  switch  in  Leg  A. 

If  C  -  ias  <  -e,  then  close  the  bottom  switch  in  Leg  A. 

If  within  the  hysteresis  band,  -e  <  C  -  iM  <  e,  then  remain  at  the  current 
state. 

This  strategy  is  most  easily  implemented  using  analog  circuitry  since  it  requires  simply  a 
fast-acting  comparator  with  hysteresis.  Problems  arise  when  high-frequency  harmonics 
corrupt  the  a-phase  current  measurements,  resulting  in  sudden  bursts  of  rapid  switching  in 
a  given  leg.  This  problem  can  be  mitigated  by  either  increasing  the  hysteresis  band  or  by 
filtering  the  current  measurement.  Both  solutions  result  in  steady-state  error  between  the 
commanded  and  actual  quantities,  where  the  level  of  error  is  oftentimes  a  complex  function 
of  load.  Additional  problems  result  when  considering  the  implementation  with  the  PUC- 
GENO.  For  hysteresis  control  for  a  given  output  sample,  the  inverter  leg  will  either  be 
"full  on"  or  "full  off."  Thus,  the  minimum  width  pulse  derivable  from  a  given  inverter  leg 
is  fixed  by  the  sampling  period.  Therefore  a  very  high  sampling  frequency  is  required  in 
order  for  the  inverter  to  respond  fast  enough  to  track  the  desired  waveform  with  sufficient 
accuracy  and  without  introducing  considerable  harmonics.  As  a  result,  bang-bang  PWM  is 
more  appropriately  implemented  using  analog  circuitry  where  the  sampling  rate  of  a  DSP 
board  is  not  at  issue.  Owing  to  this  inherent  limitation,  hysteresis  modulation  was  not 
considered  any  further  in  this  effort. 


2.  Current  Control  Algorithms 

A  PWM  Voltage-Source-Inverter  (PWM-VSI)  can  be  made  to  appear  as  an  adjustable 
current  source  by  having  the  modulating  signals  derived  from  a  fast-acting  current  control 
loop.  The  3-phase  output  voltages  are  specified  within  the  admissible  limits  described  in 
the  previous  section  in  order  to  minimize  the  error  between  the  desired  and  actual  currents. 
This  report  will  consider  two  current  control  topologies  and  underscore  the  advantages  and 
disadvantages  of  each.  In  addition  to  outlining  the  algorithm  required  for  each,  analysis 
tools  useful  for  establishing  feedback  gains  and  system  performance  are  derived  and 
illustrated  using  a  basic  example.  Both  algorithms  are  based  on  reference  frame  analysis. 
Details  of  such  analysis  is  documented  in  some  degree  in  this  report.  The  reader  is 
referred  to  supplementary  resources  for  additional  information. 


a.  Implementation  in  the  Stationary  Reference  Frame 
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The  Stationary  Reference  Frame  (STA-RF)  curent  control  is  illustrated  in  Figure  4 
where  K;  and  Kp  represent  integral  and  proportional  gains.  The  signals  iqs*  and  idf  are  the 
desired  qd-currents  in  the  STA-RF.  The  currents  iqss  and  idss  are  the  actual  qd-currents  in  the 
STA-RF  and  are  related  to  the  actual  abc-load  currents  by 


iqs  'as 

,s  J_  •  _2_  • 

**  “  73  las  73  lbs 


(33) 

(34) 


where  the  fact  that  the  phase  currents  must  sum  to  zero  in  a  three-wire  wye-connected  load 
has  been  used  to  simplify  (33)  and  (34)  so  that  only  two  current  measurements  are 
required.  For  a  delta-connected  load,  the  ab-quantities  above  would  be  the  line  currents 
which  must  also  sum  to  zero.  Expressions  for  the  signals  Vrq  and  Vrd  are  derived  from 
Figure  4  as 


V4  =  Kp  ( iqr  -  i£)  +  K,J  (iqf  -  iqss)  dt  (35) 

Vr^  =  Kp  ( U*  -  i* )  +  Ki  J  (idss*  -  idss)  dt  (36) 

These  signals  are  then  either  transformed  out  of  the  STA-RF  to  realize  the  balanced  set  of 
modulating  signals  which  would  be  the  inputs  to  STPWM  or  passed  directly  to  the  SVM 
algorithm  as  the  desired  space  vector.  For  STPWM, 


Vmod,a 
V  mod,b 

Vmod,c 


Vs 

vrq 


_  JL  ys 
2  Vrq 


=  -V, 


mod, a 


V3 
2 

Wnod,b 


rd 


(37) 

(38) 

(39) 


If  the  phase  angle  of  the  currents  is  not  critical  and  simply  a  commanded  amplitude,  I^p, 
and  frequency,  coe,  are  required,  then  iqss*  and  id^*  may  be  derived  from 


iqs  —  tamp  COS0e 


(40) 


ids  —  tamp  sin0e 


(41) 


where 


0e  =  <»e  t 


(42) 


and  ranges  from  zero  to  2k  radians. 

The  STA-RF  control  is  straightforward  to  implement  but  has  a  basic  limitation  which 
is  exposed  upon  analysis.  We  will  next  set  forth  the  analysis  of  this  architecture  then 
consider  the  selection  of  control  gains  for  the  case  of  a  balanced  3-phase  RL-load.  We 
will  employ  a  vector  notation  where  possible  to  condense  the  development.  Notationally, 
the  modulating  signals  are  related  to  the  control  signals  by 
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in  the  over-modulation  range.  The  STA-RF  phase  voltages  are  related  to  ^abcs  by 

VqdOs  =  K,*?**  (48) 

Therefore  if  we  substitute  (43)  into  (45)  and  the  result  into  (48),  we  find  that 

Vqsd0s  =  G  Vrqd0  (49) 

The  entire  system  is  most  conveniently  analyzed  in  the  Synchronous  Reference  Frame 
(SYN-RF)  since  the  forcing  functions  are  constants  in  that  frame.  As  a  result,  the  steady- 
state  variables  are  constants  and  steady-state  analysis  is  simplified.  Upon  introducing  state 
variables  Xqs  and  xd\  we  can  write  the  modeling  equations  for  the  STA-RF  implementation 


as 

pXq  =  (iqs*  -  iqs)  (50) 

pxds  =  or  -  u)  (5i) 

VqSs  =  G  K;  xqs  +  G  Kp  (iqs*  -  iqSs)  (52) 

Vdss  =  G  Ki  xds  +  G  Kp  (ir  -  i(k)  (53) 

The  equations  are  transformed  into  the  SYN-RF  by  considering  the  mapping 

Vqed0s  =  s  Ke  Yqd0s  (54) 

where 

cos0e  -sin0e  0 

sKe  =  sin0e  cos  0e  0  (55)' 

0  0  1 


or  we  can  just  consider  the  first  two  rows  and  columns  since  the  zero-sequence  quantities 
are  all  zero  in  this  formulation.  Equations  (50)-(51)  are  rewritten  in  vector  form  as 

p^s  =  (?r  -  t,dS)  (56) 

If  we  substitute  the  transformation  between  reference  frames,  (56)  becomes 

P  [(SKe)_1^qds]  =  (SK®)-1  (Tqeds*  -  T^)  (57) 

Expanding  using  the  product  rule, 

OK6)-1  p^s  +  p  (sKer^s  =  (sKe)-1  (Tqr  -  ?qi)  (58) 

and  simplifying  by  multiplying  through  by  sKe 

Ptq%  =  -  sKe  p  (sKe)-^eds  +  (tqed;  -  Tqeds)  (59) 

Upon  performing  the  matrix  multiplication  and  expanding,  we  find  that 

pXqe  =  -COeXde  +  (iqs*  -  iqs) 


(60) 
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p  X,f  =  <DeXq  +  (ids*  -  ids) 

(61) 

The  voltage  equations  (52)-(53)  are  similarly  vectorized 

=  G  K;t,ds  +  G  Kp  (Tqsd;  -  Tqsds) 

(62) 

and  transformed  into  the  SYN-RF 

(sKe)-1  fqfc  =  G  Kj  (sKe)_I  +  G  Kp  (sKer1  (Tq^‘  -  ?qds) 

(63) 

to  yield 

7^3  =  G  Ki^qds  +  G  Kp  (fq5s*  -  tqds) 

(64) 

which  expand  to 

vqes  =  G  K;  xqe  +  G  Kp  (iqf  -  iqes) 

(65) 

vdes  =  G  K;  xde  +  G  Kp  ( i<k*  -  ids) 

(66) 

Equations  (59), (60), (65)  and  (66)  represent  the  STA-RF  control  in  the  SYN-RF.  The  next 

step  is  to  formulate  the  equations  of  the  proposed  load  in  the  SYN-RF. 
is  described  by 

A  3-phase  RL-load 

Vabcs  =  rs  *abcs  ^  P  *abcs 

(67) 

where  rs  is  the  phase  resistance  and  Ls  is  the  phase  inductance, 
equations  into  the  SYN-RF  yields 

Transforming  these 

piqes  =  -  ^  iqs  -  Weil  +  Vqe$ 

(68) 

P  i<§  =  -7-  ids  +  iqs  +  T-  W&  .  (69) 

i-'S 

Therefore,  combining  (59),(60),(65),(66),(68)  and  (69)  results  in  the  following  state  space 
system  description  in  normal  form: 


which  is  in  the  form 

pj?  =  Ast  +  BsTt 

where  As  is  the  system  matrix,  t  is  the  state  vector  and  t  is  the  vector  of  inputs.  Equation 
(70)  may  be  used  to  determine  system  eigenvalues  for  various  choices  of  Kp  and  K,  and 
may  also  be  used  to  determine  the  steady-state  error.  The  steady-state  error  is  found  by 
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setting  the  derivative  in  (71)  equal  to  zero  and  solving  for  the  state  vector  t. 

t  =  -  ( As )- 1  Bs  i?  (72) 

The  difference  between  the  derived  operating  points  in  (72)  and  the  values  iqs*  and  ides* 
provides  a  measure  of  the  accuracy  of  the  current  control.  To  illustrate  the  design  process, 
consider  the  sample  parameter  values  listed  in  Table  9.  The  governing  system  eigenvalues 
may  be  evaluated  for  a  variety  of  gains.  In  addition,  suppose  a  desired  amplitude  of  10A 
is  specified  for  the  current  (iqf=lOA  and  idf=0A)  then  the  steady-state  error  may  be 
evaluated  as  well.  Sample  results  are  presented  in  Table  10. 

Table  9.  Parameters  for  STA-RF  Current  Control  Example 


rs 

Ls 

vdc 

At 

25  a 

20  mH 

850  V 

10  V 

377  rad/sec 

Table  10.  Eigenvalue  Analysis  Results 


Kp 

Ki 

eigs  1&2 

eigs  3&4 

pk  amp. 

5 

100 

-18+/-j377 

-11857+/-j377 

8. 96  A 

5 

1000 

-181+/-j377 

— 1 1693-+7— j377 

9.26A 

1 

250 

-166+/-j377 

-3209+/-j377 

1.22k 

0.1 

250 

-612+1- fill 

-790+/-j377 

7.96A 

0.01 

500 

-636 +/-j  1188 

-636+/-j434 

10.24  A 

0.01 

1500 

-636+/-j2045 

-636+/-jl291 

10.33A 

The  final  values  of  gains  in  the  table  were  then  used  in  a  detailed  ACSL  simulation  of  the 
system.  The  STPWM  algorithm  was  simulated  in  detailed.  Initially,  the  q  and  d-stationary 
reference  frame  desired  currents  were  set  to  an  amplitude  of  5A  then  around  0.06sec  the 
amplitude  was  stepped  up  to  10A.  The  results  are  depicted  in  Figure  5.  In  this  case,  the 
limitations  of  the  STA-RF  controller  are  not  apparent;  however,  what  is  obvious  is  that  the 
accuracy  is  parameter  dependent  and  as  the  load  changes,  the  steady-state  error  will  vary. 
In  addition,  the  above  table  summarizes  amplitude  information  and  does  not  address  phase 
error  which  is  critical  in  a  vector-controlled  induction  machine  drive.  The  steady-state 
error  issue  results  from  the  fact  that  we  are  attempting  to  regulate  sinusoidal  quantities.  In 
the  next  section,  the  control  action  is  moved  to  a  reference  frame  where  the  desired 
quantities  are  contants.  The  result  is  that  the  integral  action  in  the  Proportional-Integral 
(PI)  controller  facilitates  achieving  zero  steady-state  error.  One  final  point  to  be  made  is 
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that  the  eigenvalues  and  error  listed  in  Table  10  are  not  the  only  issues  driving  the 
selection  of  the  gains.  One  must  also  be  aware  that  too  large  of  gains  will  result  in 
saturating  the  control,  driving  the  STPWM  algorithm  into  either  over-modulation  or  six- 
step  modulation.  The  results  of  employing  too  large  of  control  gains  are  illustrated  in 
Figure  6  where  immediately  following  the  step  increase  in  desired  current,  the  actual 
current  tracking  is  unacceptable. 


Figure  5.  STA-RF  Current  Control  Waveforms 
(bot:  commanded  q-axis  current,  top:  actual  q-axis  current) 


Figure  6.  STA-RF  Current  Control  Waveforms  (overmodulation) 
(bot:  commanded  q-axis  current,  top:  actual  q-axis  current) 
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b.  Implementation  in  the  Synchronous  Reference  Frame 

The  Synchronous  Reference  Frame  (SYN-RF)  curent  control  is  illustrated  in  Figure  7 
where  Kiq  and  Kid  represent  the  integral  control  gains  and  Kp„  and  Kpd  represent  the 
proportional  control  gains.  The  signals  iqs  and  ids  are  the  desired  qd-currents  in  the  SYN- 
RF.  The  currents  iqes  and  i<£  are  the  actual  qd-currents  in  the  SYN-RF  and  are  related  to  the 
actual  abc-load  currents  by 


IqdOs  =  K'tbcs 


(73) 


where 


v-i 

s  3 


cos  0e  cos(  0e  -  )  COS(  0e  +  —■  ) 

2  71  2  7C 

sin  0e  sin(  0e - —  )  sin(  0e  +  ) 


2 


2 


(74) 


The  modeling  equations  for  the  SYN-RF  control  are  derived  directly  from  Figure  7  as 


P  Xq6  —  ( iq^  iqs ) 
pxd  =  (idf  -  ids) 

Vqes  =  G  Kiq  Xqe  +  G  Kpq  (iqf  -  ) 

v£  =  GKidxde  +  G  Kpd  ( i<£*  -  i«£) 


(75) 

(76) 

(77) 

(78) 


V, 


where  G  is  the  gain  of  the  PWM  method  (G=~  for  STPWM  in  the  linear  range). 

Incorporating  these  equations  together  with  (68)  and  (69)  results  in  the  following  state 
space  normal  depiction  of  the  SYN-RF  control. 


Xd 


-rs  -  G Kpq 

-®e 

GKiq 

0 

gJEei 

0 

Ls 

,  4 

1 

4 

-rs  -  GKpd 

0 

GKid 

iqs 

0 

G^- 

4 

Ls 

ids 

Ls 

-1 

0 

0 

0 

xqe 

+ 

1 

0 

0 

-1 

0 

0 

xd 

0 

1 

xqs 

•  e51 
Ids 


(79) 


which  is  in  the  form 

pt  =  Aelt  +  Be  l? 


(BO) 


One  outright  advantage  of  the  SYN-RF  implementation  is  that  in  the  steady  state  (p*  =  0), 
there  is  no  error;  that  is,  iqes  =  iqf  and  i£  =  i£".  The  "cost"  of  the  implementation  is  that 
angle-dependent  transformations  are  required  to  map  the  SYN-RF  control  signals  back  into 
abc-modulating  signals  and  the  actual  sensed  abc-currents  into  the  required  qd-feedback 


currents. 


-29- 


.  e 


Figure  7.  Synchronous  Reference  Frame  Current  Control  Implementation 

A  simple  design  illustration  is  presented  next.  Considering  the  parameters  listed  in 
Table  9,  eigenvalues  may  be  computed  for  a  variety  of  control  gains.  A  procedure  for 
placing  the  eigenvalues  for  a  related  problem  is  described  in  the  next  section  and  can  be 
readily  extended  to  this  case  study.  The  locations  of  the  eigenvalues  (or  poles)  in  the 
complex  plane  directly  influences  transient  response.  The  real-part  of  the  dominant 
eigenvalue  (closest  to  the  j  co  axis)  is  closely  related  to  the  settling  time  while  the  ratio  of 
the  imaginary  to  the  real-part  dictates  the  overshoot.  Transfer  function  zeros  in  the 
neighborhood  of  these  dominant  eigenvalues  will  also  impact  the  transient  response  and 
must  be  considered. 

Since  the  integral  action  guarantees  zero  steady-state  error  for  a  constant  input,  we  do 
not  need  to  consider  the  impact  of  the  gains  on  the  resultant  steady-state  current  amplitude. 
Table  11  documents  how  the  eigenvalues  transition  in  the  complex  plane  for  different 
gains.  Once  again,  an  effort  should  be  made  to  minimize  the  gains  in  order  to  minimize 
the  control  effort  and  avoid  overmodulation.  The  current  control  design  was  validated  in 
ACSL  using  a  detailed  representation  of  the  STPWM.  Simulation  results  are  presented  in 
Figures  8  and  9.  Initially,  the  load  operates  with  iqes*=5A  and  ides*=0A.  As  the  curves 


indicate,  there  is  no  steady-state  error.  At  about  0.06sec,  iqes*  is  stepped  to  10 A.  The 
output  current  responds  quickly  as  predicted  by  the  eigenvalues,  distortion  is  avoided  as 
overmodulation  is  not  entered  due  to  the  suitably  small  selection  of  Kp,  and  there  is  zero 
steady-state  error.  Despite  being  somewhat  more  involved  to  implement,  the  superior 
performance  and  robustness  of  the  SYN-RF  current  control  make  it  more  attractive  than  its 
STA-RF  counterpart. 


Table  11.  SYN-RF  Design  Example 


KP 

Ki 

Eigs  1&2 

Eigs  3&4 

0.05 

1000 

-580 +/-j  1119 

-776 +/-j  1496 

0.01 

1000 

-545 +/-j  1140 

—726  +/ — jl517 

0.01 

2000 

-575 +/-j  1783 

—696  +/— j2 1 60 

0.01 

250 

-384+/-j287 

-888+/-j664 

0.01 

100 

-166+/-j66 

— 1 106  +/ — j443 

3.  Voltage  Control  Algorithms 

On  board  a  naval  ship,  it  may  be  more  common  to  require  an  ARCP  inverter  to 
establish  and  regulate  a  balanced  set  of  three-phase  voltages.  Since  it  is  generally  required 
that  these  be  sinusoidal  voltages  with  very  little  harmonic  distortion,  the  normal  PWM 
output  voltages  are  processed  through  a  passive  filter  before  being  applied  to  a  load.  In 
order  to  illustrate  a  high-performance  scheme  involving  both  feedforward  and  feedback 
components,  consider  the  diagram  illustrated  in  Figure  10.  Note  that  the  load  is  once  again 
assumed  to  consist  of  a  resistor  and  an  inductor.  The  inverter  is  coupled  to  the  load 
through  an  LC-filter.  The  governing  equations  for  this  network  in  terms  of  system 
variables  are  given  in  vector  form  as 


—  YaBCii  Lf  P  lLf,abc  ^L,abc  0 

.  T>  __  n 

Cf  p  vLabc  =  lLf,abc  ”  ^L.abc 

-  ^L,abc  +  RL?L,abc  +  LL  P  lL,abc  =  0 


(81) 

(82) 

(83) 


where 


YaBCii 


VAn 

VBn 

VCn 


(84) 
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(86) 

(87) 


Figure  10.  Inverter  Voltage  Control  Output  Filter  and  Load 


It  has  been  shown  in  previous  sections  that  analysis  in  the  Synchronous  Reference  Frame 
(SYN-RF)  is  advantageous  in  modeling  inverter  systems  while  considering  only  the 
fundamental  quantity.  This  holds  because  sinusoids  become  constants  in  the  SYN-RF 
implying  that  all  state  variables  are  constants  in  the  steady  state.  If  we  consider  balanced 
operation  where  all  zero-sequence  variables  are  identically  zero,  we  can  transform  (81)-(83) 
into  the  SYN-RF  and  arrive  at 


(88) 

(89) 


-  V£d  + 


(90) 


QSE 
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Figure  9.  SYN-RF  Current  Control  a-Phase  Current  Waveforms 

P  VLe,d  =  ®e  Vuq  +  -pr  iLf,d  ~  ~  (91) 

M  M 

P  il!q  =  -®e  iLe,d  +  “  vl!q  -  il!q  (92) 

P  iLed  =  iLeq  +  J-  v£d  -  iLe,d  (93) 

lL 

Thus,  our  system  inputs  are  VQn  and  Vpn  which  are  arrived  at  from  manipulating  the  inverter 
switching  via  either  STPWM  or  SVM.  In  terms  of  our  original  quantities,  v<fn  and  v^n  are 
the  SYN-RF  transformed  versions  of  v^bcd- 

We  desire  to  regulate  the  output  load  voltages  represented  by  vLeq  and  vLed.  Let’s 
assume  that  we  can  monitor  the  instantaneous  values  of  the  output  voltages  (vLeq  and  vLed), 
the  output  currents  (i^q  and  iLed),  the  filter  inductor  currents  (iLefiq  and  i^d),  and  the  input  dc 
voltage  (Vdc).  That  is,  assume  that  the  appropriate  abc-quantities  can  be  sampled  and  then 
transformed  into  the  SYN-RF  at  the  front-end  of  the  DSP  control  implementation. 
Therefore,  assuming  balanced  operation,  the  proposed  scheme  requires  seven  variables  to 
be  sampled. 
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Before  we  delve  into  the  control  architecture,  consider  that  the  desired  output  voltages 
assume  the  following  form 


(94) 


* 

V0  COS(COet) 

vL,a 

vL*,b 

= 

,,  .  2k  , 

V0  cos(  toe  t  -  — ) 

vL,c 

,  2k  . 

V0  cos(coet+  —  ) 

Transforming  this  balanced  set  into  the  SYN-RF  yields 


VL,q 

VL,d 

II 

0 

vl,o 

.0  . 

(95) 


We  can  therefore  view  our  setpoint  as  being  vLe*.  Also  since  the  load  voltage  is  across  the 
filter  capacitor,  it  follows  that  the  desired  capacitor  current  is  given  by 

ic,abc  =  Q  P  ^L.abc 

or  transforming  into  the  SYN-RF 


lc,qdO 


=  Qp^i 


L,qdO 


©e  Q  vL,d 

-we  cf  vLe; 

0 


(97) 


The  proposed  control  architecture  is  illustrated  in  Figure  11.  We  will  first  describe 
the  flow  qualitatively  in  order  to  provide  the  reader  with  the  basic  idea.  First,  there  is  an 
inner  fast-acting  current  control  loop  surrounded  by  a  slower-acting  voltage  control  loop. 
The  voltage  control  loop  consists  of  Proportional-Integral  (PI)  regulators  for  both  the  q  and 
d-axis  output  voltages.  The  outputs  of  these  regulators  provide  for  transient  changes  in  the 
commanded  filter  capacitor  currents.  The  commanded  filter  capacitor  currents  also  consist 
of  a  feedforward  term  which  programs  in  the  required  steady-state  currents.  Next,  the  load 
curent  is  added  to  the  capacitor  current  to  establish  a  commanded  value  of  inductor  filter 
current.  The  inner  current  control  processes  the  error  between  the  commanded  filter 
inductor  current  and  the  actual  filter  inductor  current  to  generate  q  and  d-feedback  derived 
signals.  The  desired  inverter  output  voltages  are  found  by  adding  the  feedback  signals  to 
feedforward  components,  dictated  by  the  commanded  capacitor  voltages  and  the  voltage 
drops  across  the  filter  inductances.  These  signals  then  are  either  transformed  to  the  STA- 
RF  for  SVM  or  to  the  abc-realm  for  STPWM. 

Certainly  there  is  nothing  trivial  about  this  formulation.  Let’s  step  back  and  walk 
through  it  one  more  time  with  the  aid  of  the  equations. 


Step  1.  Perturbation  signals  are  derived  from  the  error  in  the  commanded  and 
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actual  output  load  voltages: 

icap.q  =  Kpvq  ( VLe,‘q  -  VLeq  )  +  Kivq  J  (  VLe*q  -  VL*q  )  dt  (98) 

iceap,d  =  Kpvd  ( vLe;d  -  v£d )  +  Kivd  J  ( vLe;d  -  vLed )  dt  (99) 


Note  from  (97)  that  with  the  commanded  voltage  signals  constant,  the  steady-state 
capacitor  filter  currents  are  given  by 


Te* 


c,qdO 


t°e  Q  Vjf  d 
-COeCfVLeq 
0 


(100) 


Step  2.  Incorporate  feedforward  action  to  form  the  desired  filter  capacitor  currents: 

ic?q  =  icap,q  +  ®e  Q  V^d  (101) 

ic?d  =  icap.d  -  «eCf<  (102) 

Thus,  ic|pq  and  iceapd  "correct"  the  feedforward  components  and  are  necessary  to 
compensate  for  imprecise  knowledge  of  Q  and  for  loss  mechanisms  ignored  in  the 
analysis. 


Step  3.  Form  the  desired  filter  inductor  currents  from  the  desired  filter  capacitor 
currents  and  the  measured  load  currents  (in  the  SYN-RF): 


•  e*  _  :  e*  ,  :  e 

JLf,q  —  Jc,  q  !L,q 

iLtd  =  kj  +  iL,d 

These  signals  become  inputs  to  the  inner  current  control. 


(103) 

(104) 


Step  4.  Perturbation  signals  are  derived  from  the  error  in  the  commanded  and 
actual  filter  inductor  currents: 


v4,q  =  Kpcq  (iL?q  -  iftq)  +  Kicq  j  (iLtq  -  iLf,q)  dt  (105) 

vfbk,d  =  Kpcd  (iud  ~  iu.d)  +  kjcd  j  ( iu,d  -  iu,d )  dt  (106) 


Rewriting  (88)  and  (89),  we  have  that 

Vq„  =  V^q  +  Lf  p  ijjf,q  +  C0eLfijjr  d  (107) 

vDn  =  vL,d  +  LfpiLf  d  -  ®e  Lf  iLef  q  (108) 

Thus  in  the  steady  state  we  can  ignore  the  derivative-of-current  terms  and  generate  a 
feedforward  component  for  the  desired  line-to-neutral  voltage  across  the  filter  capacitor  and 
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Figure  11.  Voltage  Control  Block  Diagram 
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inductor.  We  accomplish  this  by  measuring  the  actual  load  voltages  and  using  the 
previously  derived  commanded  filter  currents. 

Step  5.  Form  the  commanded  line-to-neutral  voltage  across  the  filter  capacitor  and 
inductor  from  the  steady-state  feedforward  component  using  the  commanded  filter 
inductor  currents  and  the  feedback  component  generated  in  Step  4: 

VQn  =  vL,q  +  Lf  iuy  +  Vfbtq  (109) 

vDn  =  vL,d  “  °>eLfiLf,q  +  Vfl^d  (HO) 

Thus,  v4_q  and  "correct"  the  feedforward  components  to  compensate  for 
uncertainty  and  non-idealities. 

Step  6.  Use  vq*  and  vgl  to  generate  the  required  STPWM  or  SVM  signals.  That  is, 
either  transform  these  quantities  to  abc-quantities  for  use  as  modulation  signals  in 
STPWM  or  transform  them  to  STA-RF  quantities  for  application  in  SVM. 

Synthesizing  the  control  gains  in  a  systematic  manner  is  a  non-trivial  endeavor, 
though  the  process  can  be  simplified  by  making  some  basic  observations  and  assumptions. 


The  filter  inductor  current  dynamics  are  given  by 

Pil?f,q  =  -  “e  iLf.d  +  7-VQn  -  J~  VL,q  (HI) 

Lf  L{ 

piu,d  =  “e  iLf.q  +  ^  VDen  -  ^  VL,d  (H2) 

where  the  actual  applied  SYN-RF  voltages  out  of  the  inverter  are  related  to  the  control 

signals  by  the  gain  G  as  shown  in 

VQn  =  G  Vq„  =  G  (  VLe_q  +  C0e  Lf  i&*d  +  V^)  (113) 

VD„  =  GvDen*  =G(V£d  -  coeLfiLtq  +  VfgM)  (114) 

The  dynamics  introduced  by  the  current  control  integral  compensators  are  described  by 

pxceq  =  iuq  -  iu,q  (H5) 

P  xcd  =  ^Lf,d  -  il!f,d  (116) 

and 

Vfbk,q  =  kpcq  (ii^q  —  'Lf.q )  +  K;Cq  XCq  (117) 

vftk,d  =  kpcd  ( ijjF.d  -  iLf,d)  +  Kicd  Xcd  (118) 


Therefore,  if  we  choose  G  in  (1 13)-(1 14)  to  be  unity  and  substitute  the  result  into  (111)- 
(112),  the  filter  current  dynamics  become  simply 
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pil?f,q  .=  -  ®e  iLef,d  +  ®eiLf*d  +  ~  (»Uq  “  W  +  Lf  ** 


(119) 


KPcd  n  e*  _  ;e  \  ,  JE“i  v  « 


piLf,d  =  <Muq  "  “e  lUq  +  “  (»&  ~  »Ud)  +  Lf  X<^d 

Combining  (115)-(116)  and  (119)-(120)  yields  the  state  space  description: 


(120) 


Kjcd  ‘Lf,q 

Lf  iiff.d 


(121) 


Instead  of  using  trial-and-error,  the  eigenvalues  for  (121)  may  be  uniquely  specified  using 
the  four  available  gains.  To  illustrate  how  this  is  accomplished,  consider  the  characteristic 

equation 

^4  +  +  +  +  +  <ft?)xa  (122) 


4  „  thu  +  hsi.xi  +  +  +5a. hsL  +  ^)X2 


{  ^icci  ,  ^icd  ^icd-  =  q 

Lf  Lf  Lf  Lf  Lf  Lf 

Equation  (122)  can  be  compared  with  a  desired  characteristic  equation  formulated  from  the 
desired  eigenvalue  locations.  Assume  that  the  desired  equation  is  given  by 

X4  +  d3  X3  +  d2  X2  +  d,  X  +  do  =  0  (123) 

If  we  let  Xl  =  ^-,  x2  =  ^,  x3  =  ^,  and  *4  =  ^,  then  we  can  coefficient  match 
Lf  Lf  Lf  f 

between  (122)  and  (123)  and  formulate  the  following  set  of  four  nonlinear  equations  in 
four  unknowns. 


fj  =  xj  +  x2  -  d3 

f2  =  XfX2  +  X3  +  x4  +  Cl)e2  -  d2 
f3  =  XlX4  +  X2X3  -  d! 
f4  =  x3x4  -  do 


(124) 

(125) 

(126) 
(127) 


where  we  seek  to  find  values  of  ’x’  where  each  function  T  is  equal  to  zero.  We  can  solve 
this  problem  by  using  a  Newton-Raphson  approach  where  we  define  the  vector-valued 
nonlinear  vector-function 


(128) 
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our  vector  of  unknowns 


Xl 

X2 

X3 

X4 


(129) 


and  the  Jacobian  matrix 

'  1  1  o  o' 

X2o  Xl0  1  1 

F  =  X4o  X3o  X2o  Xio 

0  0  X^  X3o 

where  x]o  through  x^  represent  the  current  guesses  of  the  unknowns  X!  through  x4.  We  can 
solve  for  the  unknowns  by  iteratively  applying 


(130) 


*n+1  =  %  -  F (*„)-' ?(tn) 


(131) 


An  example  will  be  illustrated  at  the  end  of  this  section. 

The  voltage  loop  may  then  be  designed  by  assuming  that  the  current  dynamics  are 
much  faster  than  the  voltage  control  dynamics  so  that 


iu,q  =  iuq 

iud  =  iua 

The  dynamics  introduced  by  the  voltage 


(132) 

(133) 

control  integral  compensators  are  described  by 


P  xveq  =  vLe*  -  vLe>q 
px,d  =  v£*d  -  vLed 

The  current  loop  being  assumed  ideal  allows  us  to  rewrite  (90)-(91)  as 

l  *  l 

P  VLe,q  =  ^  itf,q  “  'Uq 

P  VLe,d  =  C0e  v£q  +  ~~  iLM  -  ^  iL,d 


(134) 

(135) 


(136) 

(137) 


which  can  be  reduced  upon  inspection  of  the  control  algorithm  since 

il$q  =  i£q  +  ®eCfVLed  +  Kpvq  ( <  -  VLe>q )  +  Kjvq  Xveq  (138) 

iu,6  =  iud  -  ®e  Cf  vLeq  +  KpVd  (vLe‘d  -  vLe,d)  +  Kivd  xved  (139) 

Thus,  substituting  (138)-(139)  into  (136)-(137)  then  augmenting  with  equations  (134)-(135), 
we  arrive  at  One  can  quickly  ascertain  that  the  system  matrix  is  identical  in  form  to  the 
system  matrix  in  (121).  Therefore,  the  eigenvalue  assignment  problem  is  the  same:  solve 
an  identical  set  of  four  nonlinear  equations  using  the  previously  described  Newton-Raphson 
algorithm. 


-40- 


P 


VM 

VL,d 

xvq 

xvd 


Kpvq 

Cf 

“e 

-1 

0 


—  ®e 
Kpvd 

Cf 

0 

-1 


Kjyq 

“cT 

0 

0 

0 


0 

Kjvd 

"o' 

0 

0 


Xvq 

Xvd 


CDe 


+ 


1  0 
0  1 


(140) 


Next,  let’s  illustrate  the  procedure  with  a  design  example.  Consider  the  parameters 
listed  in  Table  12. 


Table  12.  Parameters  for  Inverter  Voltage  Control  Example 


Rl 

Ll 

Q 

Lf 

25  £2 

20  mH 

0.1  mF 

10.1  mH 

Vdc 

A, 

u 

G>e 

850V 

425  V 

5  kHz 

377  rad/sec 

A  set  of  MATLAB  script  files  were  written  to  perform  the  calculations.  First,  the  four 
current  control  eigenvalues  were  specified.  For  convenience,  all  four  were  placed  at  the 
same  location,  comfortably  in  the  left-half-plane,  yet  not  too  far  to  interfere  with  the 
switching  frequency  2jt5000rad/sec. 

^-1,2,3, 4  =  “  3000 

Next,  the  current  gains  were  found  by  applying  the  Newton-Raphson  algorithm.  The 
voltage-loop  eigenvalues  were  then  specified.  We  chose  to  locate  them  a  factor  of  five 
times  closer  to  the  jo>  axis  so  that  the  current-control  dynamics  would  be  faster.  In 


particular, 


^5, 6, 7, 8  —  —  600 

The  voltage-loop  gains  were  then  found  by  applying  the  Newton-Raphson  algorithm.  Both 
sets  of  gains  are  recorded  in  Table  13.  The  eigenvalues  of  the  overall  system  were  then 
established  without  using  any  loop-interaction  approximation  yielding  the  following  results: 


-2649+/-j2699 
-2472 +/-j  1926 
—825  j415 

— 459  -t-/ — j  100 
— 653 
-1040 


Note,  the  actual  eigenvalue  locations  are  somewhat  different  than  the  desired  locations. 
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This  occurs  because  of  the  approximation  that  the  current  loop  is  infinitely  faster  than  the 
voltage  loop.  Had  we  located  the  current-control  eigenvalues  further  into  the  left-half¬ 
plane,  the  resultant  eigenvalues  would  more  closely  track  the  desired  values.  Nonetheless, 
the  above  concentration  is  approximately  where  we  want  them  and  so  we  should  expect 
acceptable  transient  response  (though  there  are  some  zeros  near  to  the  j  ©  axis  which  tend 
to  increase  the  overshoot). 


Table  13.  Multi-Loop  Control  Gains 


Kpcq 

Kjcq 

Kpcd 

Kicd 

56.79 

90900.0 

64.41 

90900.0 

Kpvq 

Kjyq 

Kpvd 

Kjvd 

0.0823 

36.0 

0.1577 

36.0 

The  next  step  in  the  procedure  is  to  validate  the  design  using  a  detailed  ACSL 
simulation.  One  of  the  previous  programs  was  modified  to  depict  the  described  algorithm. 

Initially,  it  is  assumed  that  the  filter  capacitors  are  de-energized.  As  a  result,  the 
commanded  voltage  is  ramped  up  from  OV  to  360V,  over  a  0.1  sec  interval.  That  is,  vLe*  is 
ramped  from  0V  to  360V  while  vLe*d  is  held  at  zero.  The  360V  corresponds  to  the  peak 
phase  voltage;  therefore,  we  are  ramping  up  to  a  rms  line-to-neutral  voltage  of  255V  or  an 
rms  line-to-line  voltage  of  440V.  The  0.1  sec  is  arbitrary  but  is  required  to  limit  the  current 
drawn  from  the  inverter  and  helps  maintain  the  inverter  in  the  linear  modulation  range.  At 
0.15sec,  the  load  is  step-changed  from  25Q  down  to  10Q.  The  plots  illustrated  include  the 
a-phase  commanded  and  actual  phase  voltage,  the  q  and  d-load  voltages,  the  filter  and  load 
currents,  and  the  a-phase  PWM  modulating  signal.  Note  that  the  integral  action  on  the 
qd-load  voltages  guarantees  zero  steady-state  error.  Also,  the  fast-acting  inner  current  loop 
together  with  the  feedforward  action  guarantees  rapid  response  to  the  load  change  and  only 
a  minimal  voltage  transient.  The  transient  persists  for  about  10msec  as  anticipated  with  the 
dominant  eigenvalue  located  with  a  real  part  of  approximately  -400.  The  q-load  voltage 
decreases  by  about  4V  while  the  d-load  voltage  decreases  by  about  2V  resulting  in  less 
than  a  5V  decrease  in  the  line-to-neutral  voltage,  hardly  even  noticeable  in  the  figure.  As 
is  illustrated,  the  voltages  are  very  sinusoidal  and  the  modulation  signal  remains 
predominantly  in  the  linear  modulation  range.  In  addition,  it  is  important  to  recognize  that 
the  control  design  was  performed  independent  of  the  load  parameters.  The  gains  were 
selected  based  solely  on  knowing  the  filter  parameters  and  the  desired  output  frequency. 
The  ultimate  selection  of  the  inverter  output  filter  components  will  be  dictated  by  stability 
issues  and  the  desired  level  of  harmonics. 
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The  algorithm  may  be  summarized  as  follows. 

1.  The  a  and  b-phase  load  voltages,  load  currents  and  filter  currents  are 

measured  with  appropriate  sensors. 

2.  The  q  and  d-synchronous  reference  frame  quantities  are  then  calculated. 
For  example  for  the  load  voltage,  this  may  be  shown  in  two  steps  (first 
going  to  the  stationary  reference  frame) 

VL,q  =  VL,a 

VLS,d  =  -^VL>a  -  ^VL,b 

(Then  to  the  synchronous  reference  frame) 
vLe,q  =  vLs,q  cos0e  -  v£d  sin0e 
vLe,d  =  vL,q  sin6e  +  VLS>d  COS0e 

where  0e  is  the  angle  of  the  SYN-RF  that  varies  at  coe.  This  angle  or  the 
trigonometric  functions  may  be  evaluated  by  table  lookup  as  described  by 
Floodeen. 

3.  The  commanded  load  voltage  components  are  initialized  wherein  a  ramp 
algorithm  may  be  used  to  initiate  startup. 

4.  Equations  (98)-(99)  and  (101)-(110)  are  programmed  using  the 
trapezoidal  approximation  for  the  integration  function  as  described  by  both 
Hanson  and  Floodeen. 

5.  The  appropriate  duty  cycle  signals  are  formulated  using  either  STPWM 
or  SVM. 

Thus  far,  current  control  and  a  voltage  control  algorithm  has  been  considered  for  the 
ARCP  inverter  outer  control.  The  final  element  of  this  report  will  document  an  induction 
machine  start-up  and  speed  control  algorithm  that  will  build  on  the  previous  presentation. 
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Figure  12.  Voltage  Waveforms  for  the  Voltage  Control  Example 
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Figure  13.  Current  Waveforms  for  the  Voltage  Control  Example 
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Figure  14. 
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B.  Induction  Machine  Startup  Control 

On  board  ship,  each  electrical  zone  will  have  numerous  electric  machine  loads  that  in 
turn  power  compressors,  blowers,  fans  and  other  actuators.  In  DC  ZEDS,  most  of  these 
machines  are  induction  motors  and  all  must  be  energized  by  SSIMs,  providing  variable- 
voltage  variable-frequency  control  opportunites.  It  is  well  understood  that  Direct-On-Line 
(DOL)  starting  of  an  induction  machine  typically  results  in  a  large  in-rush  current.  One 
approach  to  limiting  the  starting  current  is  to  reduce  the  applied  voltages  at  rated  frquency. 
This  has  the  disadvantage  of  also  reducing  the  available  starting  torque,  which  may  or  may 
not  be  an  issue  for  the  load.  By  varying  both  the  voltage  amplitude  and  the  applied 
frequency,  the  machine  stator  current  can  be  limited  without  compromising  the  starting 

torque. 

In  this  section,  we  wish  to  propose  two  approaches  to  implementing  an  induction 
machine  startup  control  that  allow  the  user  to  ramp  up  the  speed  of  the  machine  to  a 
desired  level  and  have  that  speed  regulated  in  the  presence  of  any  load  disturbances.  In 
addition,  it  is  desired  to  automatically  limit  the  current  impressed  through  the  SSIM 
semiconductor  switches  and  avoid  the  large  surge  of  current  from  the  dc  bus  when  the 
machine  is  starting. 

Both  approaches  proffered  use  a  multi-loop  design  where  a  slow  speed-control  loop  is 
wrapped  around  a  fast-acting  current-control  loop.  Since  the  machine  torque  is  intimately 
related  to  the  stator  currents,  we  may  view  the  outer  loop  as  generating  a  commanded 
torque  signal  and  the  inner  loop  as  programming  the  currents  necessary  to  realize  that 
desired  torque.  The  current  control  operates  in  a  manner  analogous  to  the  approach 
outlined  in  a  previous  section.  The  outer  speed  control  loop  develops  commanded  values 
of  frequency  and  current  amplitude  to  indirectly  control  the  torque  and  regulate  the 
magnetic  flux.  If  the  machine  flux  is  not  held  below  the  saturation  level  of  the  iron, 
additional  losses  are  introduced,  harmonic  distortion  is  increased,  and  airborne  noise  and 
vibration  may  become  unacceptable. 

1.  Constant  Air-Gap  Flwc  Control 

The  standard  per-phase  steady-state  equivalent  circuit  of  the  induction  machine  is 
illustrated  in  Figure  15.  The  parameters  rs  and  t'  are  the  resistances  of  the  stator  winding 
and  the  referred  rotor  winding  (referred  to  have  the  same  number  of  turns  as  the  stator). 
The  reactances  Xls  and  are  associated  with  the  leakage  flux  produced  by  the  stator  and 
referred  rotor  windings.  The  reactance  Xm  accounts  for  the  production  of  the  flux  which 
couples  both  the  stator  and  the  rotor,  the  air-gap  flux.  The  resistance  Rc  is  oftentimes 
inserted  into  the  circuit  to  account  for  core  losses  arising  from  hysteresis  and  eddy 
currents.  The  current  flowing  down  through  the  center  branch  is  termed  the  magnetizing 
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Figure  15.  Per-Phase  Steady-State  Equivalent  Induction  Machine  Circuit 


current  and  is  directly  related  to  the  air-gap  flux.  If  we  assume  that  our  fast-acting  current 
control  resembles  a  current  source  (the  electrical  dynamics  are  much  faster  than  the 
mechanical  dynamics),  then  the  stator  phase  current  is  fixed  and  neither  rs  or  Xb  will 
effect  its  value.  How  that  current  divides  into  magnetizing  current  and  rotor-referred 
current  is  dictated  by  the  machine  slip 


s  = 


(141) 


where  <oe  is  the  fundamental  frequency  of  the  applied  1^,  av  is  the  actual  rotor  speed  in 
radians/sec,  and  tor  is  the  rotor  electrical  speed  which  is  related  to  com  by  the  number  of 

p 

pole  pairs, 

For  small  values  of  slip,  the  rotor  branch  impedance  is  large  and  thus  most  of  the 
stator  current  proceeds  down  the  magnetizing  branch.  Therefore,  since  rated  magnetizing 
current  is  typically  much  smaller  than  rated  stator  current,  a  moderate  amplitude  of  1^  has 
the  potential  of  saturating  the  machine  iron.  As  a  result,  we  must  ensure  that  for  a  given 
operating  speed  and  applied  frequency  that  the  appropriate  amplitude  of  current  is  applied. 
The  magnetizing  current  is  related  to  the  stator  current  by 
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where 


rr  coe 

+  iTrXt 
_s _ ©b 

Tr  0).  , 

+ 

s  (0b 


Xrr  =  Xm  +  Xlr 


ias 


(142) 


(143) 


If  we  constrain  operation  to  occur  at  rated  magnetizing  current,  Im  =  lmja„  then 
magOUra,)2  = 


rrr  2  r«bx'v2 

(T>  + 


<->2  +  ( —  X^)2 


magCIas)2 


(144) 


©b 


dividing  each  term  by  ( —  )2  and  defining 

COb 


co„ 


©sin  =  S' 


coe-cor  coe 


ooe-cor 


COb  ©e  ©b  ©b 

where  oosin  is  the  normalized  slip  frequency  yields 


©si 

©b 


(145) 


magtVrat )  — 


(  — )2  +  (Xt)2 

^sln 


-  mag(  las  )2 


(146) 


(— )2  +  (X^)2 

®sln 


(147) 


which  can  be  rearranged  to  solve  for  the  required  stator  current  amplitude 

/  ( Xn- cosi )z  +  (rrcotr  _  ,~T  N 

ma8<1“)  “  "V  (X>,)2  T  (rX)2  ma8(I~“) 

Therefore  given  the  rated  magnetizing  current  of  the  machine,  the  applied  electrical 
frequency  coe  and  the  rotor  electrical  angular  velocity  cOf,  we  can  establish  the  required 
stator  current.  This  relationship  will  figure  predominately  in  the  proposed  control  strategy. 
Further  it  can  be  shown  that  while  operating  at  rated  flux,  the  developed  torque  is  directly 
related  to  the  slip  frequency.  In  particular, 


T„  =  3 


P  xm2 


rr  s  coen 


2  ©b  rr  2  +  ( Xjf  s  coen  y 


mag  dm,  rat ) 


where 


_©e_ 

©b 


For  scoen  small,  (148)  reduces  to 


T„  =  3  4* 


P  X^  sco, 


2  ©b  rr 


—  mag(  Imrat  )2 


(148) 


(149) 


(150) 
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But  since  scoen  =  —  — ,  it  follows  then  that  the  torque  is  approximately  proportional  to 

©b 

the  slip  frequency  we-cor  for  values  up  to  rated  torque  when  the  air-gap  flux  is  held 
constant. 

Thus,  we  have  a  mechanism  for  achieving  speed  control.  The  outer  speed  control 
loop  develops  a  commanded  slip  frequency  signal  which  directly  controls  the  torque.  By 
sensing  the  actual  rotor  speed,  this  signal  can  be  added  to  the  slip  frequency  signal  to 
derive  the  required  fundamental  frequency  of  the  inverter  modulating  signals.  In  addition, 
the  commanded  slip  frequency  can  be  used  together  with  (147)  to  specify  the  amplitude  of 
the  stator  current  required  to  maintain  the  air-gap  flux  at  its  rated  level. 

Consider  the  block  diagram  illustrated  in  Figure  16.  A  commanded  slip  frequency 
signal  is  established  by  the  outer  speed  control  loop.  The  rotor  electrical  angular  velocity 
is  added  to  this  signal  to  derive  the  desired  coe*.  In  addition,  co‘,  is  then  used  to  establish  the 
proper  stator  phase  current  amplitude  to  guarantee  rated  air-gap  flux.  This  peak  amplitude 
is  then  allocated  into  iqf  and  ig*.  The  inner  current  control  loop  then  forces  the  actual 
current  to  track  these  desired  values.  The  maximum  torque,  at  rated  air-gap  flux,  occurs  at 

cosi,max  =  ®b  rb-  (151) 

Therefore  we  can  maintain  the  torque  on  the  negative-sloped  portion  of  the  torque  versus 
speed  characteristic  by  imposing  this  limit  on  the  commanded  slip  frequency  signal.  In 
general,  we  would  probably  specify  a  lower  value  so  as  to  limit  the  amplitude  of  the 
commanded  currents  out  of  the  inverter. 

Let’s  consider  how  a  dynamic  plays  out.  Consider  a  step  increase  in  the  commanded 
rotor  speed.  The  output  signal  of  the  speed  controller,  cos*,  saturates  corresponding  to  a 
demand  of  maximum  torque.  This  is  achieved  by  commanding  the  appropriate  maximum 
value  of  current.  As  the  current  control  quickly  establishes  the  required  current,  the 
developed  torque  accelerates  the  rotor  and  the  speed  error  decreases.  Eventually,  to*,  comes 
off  the  limit,  the  demanded  stator  current  amplitude  is  reduced,  and  the  developed  torque 
decreases  until  steady  state  is  achieved. 

The  control  requires  that  two  stator  currents  be  sensed  for  the  inner  current  control 
and  that  a  rotor  speed  signal  be  derived  for  the  outer  loop  calculations.  The  bounds  on  the 
commanded  slip  frequency  and  the  associated  commanded  current  can  be  adjusted  to 
appropriately  limit  the  current  being  demanded  out  of  the  inverter. 

The  above  control  was  implemented  for  a  4-pole,  j  horsepower  squirrel-cage 

induction  machine  in  the  laboratory  at  NPS.  The  machine  has  a  rated  rms  phase  voltage  of 
120V,  rated  rms  phase  current  of  1.2A,  and  a  rated  speed  of  1670rpm.  The  parameters  of 
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current  feedback 


Figure  16.  Constant  Air-Gap  Flux  Induction  Machine 
Speed  Control  Block  Diagram 
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the  machine  as  determined  by  resistance  measurements,  no-load  and  blocked-rotor  tests  are 
summarized  in  the  following  table. 

Table  14.  NPS  Induction  Machine  Parameters 


rs  =  11.78  Q 

Xk=  10.84  a 

xm 

=  149.0Q 

rr  =  8.83  £2 

Xk  =  10.84  Q 

Rc 

=  2097  Q 

P  =  4 

J  =  0.0024  kg/m2 

B 

=  0.001  N-m-sec 

In  order  to  understand  how  the  various  control  gains  can  be  analytically  determined,  we 
need  to  set  forth  the  modeling  equations  for  the  induction  machine  in  the  synchronous 
reference  frame. 


•^iqes  =  -Cliqs  -  ©e  ids  +  ^Vqr  "  °3  “Wd?  +  C4Vqs 
ides  =  ©eiqs  -  Cjij®  +  C3©rV]/q?  +  C2^e  +  C4Vde$ 

-^Vqr6  =  C6iqes  -  C5Vq?  -  (  ©e  ~  ©r  )  V* 

4-Vdr6  =  C6i|  +  (©e~©r)Vqre  ~ 

at 

4-(0r  =  -c7©r  +  c8iqesy*  -  Cgi^i|/q're  -  c9Tl 
at 

where  the  parameter-dependent  constants  c,  through  c9  are  given  by 


c2  - 


C3 


C4  = 


C5  = 


C6  = 


(152) 

(153) 

(154) 

(155) 

(156) 


©b 

D 

IT  +  H 

(157) 

rr  ©b  Xm 

ox: 

(158) 

xm 

D 

(159) 

xX 

D 

(160) 

rr  ®b 

X„ 

(161) 

rr  ©b  Xm 

X„ 

(162) 

B_ 

(163) 
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3P2Xm  (164) 

8JcObX^ 


c9  = 


_P_ 

2J 


(165) 


The  interconnection  of  the  synchronous  reference  frame  current  control  and  STPWM  to  the 
induction  machine  representation  is  accomplished  in  a  manner  similar  to  that  described 
previously.  In  particular,  the  integral-action  in  the  current  compensators  introduce  the 
following  dynamic  equations 


(166) 


idf  -  id* 


(167) 


and  the  induction  machine  input  voltages  are  then  given  by 


vqs  =  Kpwm  Kpq  ( iqs  -iqs)  +  Kpwm  Kiq  Xq 
V(js  =  Kpwm  Kpd  (ids  —  ids  )  +  Kpwm  Ky  Xd 


(168) 

(169) 


where  KpWm  is  the  gain  of  the  PWM  modulator  which,  in  the  linear  range,  is  given  by 


is  (170) 

J^pwm  2  A, 

Here  we  have  once  again  assumed,  for  analysis  purposes,  that  the  fundamental  inverter 
output  voltages  dominate  in  writing  (168)  and  (169).  The  outer  speed  control  loop 
contributes  the  following  dynamic  equation 


(171) 


and  the  output  equation 

00s*!  =  Kps  (CO*  -  CDr)  +  KisXs  (172) 

The  gain  of  the  block  implementing  the  relationship  between  the  commanded  slip 
frequency  and  the  commanded  stator  current  amplitude  may  be  approximated  by  the 

following, 


K, 


las, max  Im,rat 


(173) 


Before  developing  expressions  used  to  determine  feedback  gains,  let  s  consider  how 
Iaspk  may  be  employed  by  the  inner  current  control.  Recall,  the  inner  current  control 
utilizes  the  inputs  iqf  and  idf.  We  could  simply  set  one  equal  to  zero  and  set  the  other 
equal  to  pk.  The  dynamic  that  results,  however,  is  unacceptable.  For  instance,  if  we  set 
i<g*  equal  to  zero,  then  iqes*  would  be  responsible  for  both  maintaining  constant  air-gap  flux 


-  53  - 


and  adjusting  the  electromagnetic  torque.  A  better  approach,  somewhat  mimicking  what  is 
done  in  vector  control,  would  be  to  use  i£*  to  establish  a  baseline  value  of  air-gap  flux  and 
use  iqf  to  control  the  torque.  This  may  be  accomplished  by  setting 


<ds  =  I, 


■m,rat,pk 

iqf  =  Vlr 


;,pk  *ds 


(174) 

(175) 


The  modeling  equations  described  above  are  nonlinear.  In  order  to  investigate  operating 
point  stability,  we  must  linearize  the  equations.  Upon  noting  that, 


Acos*  =  Kps(Aco*  -  Acor)  +  KisAxs  (176) 

where  A  indicates  a  perturbation  quantity  away  from  an  assigned  operating  point.  It 
follows  that 


Aooe  =  Aco*!  +  Acor 


(177) 


Also,  noting  that  we  may  approximate  our  nonlinear  relationship  between  the  commanded 
slip  and  desired  stator  current  amplitude  as 


Aas,pk 


__  * 
=  K,cosl 


(178) 


Therefore, 

Alaspk  =  Kj  A(0*|  (179) 

Finally,  since  our  expression  for  iqes*  is  nonlinear,  linearizing  it  yields 


Aiqs 


*as,pko 


Vlas.pko  idso 


2  AIas>pk 


=  K, 


(180) 


where  I  k0  is  the  peak  value  of  the  stator  current  at  the  selected  operating  point. 
Therefore,  linearizing  the  machine  equations  and  substituting  the  above  results,  we  arrive  at 
the  following  model: 
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A£iml  =  [  Aiq'  Ai|  Avq'rc  A\|rd're  A(0r  Axqe  Ax/  Axs] 
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Biml 


-itg0Kps+d3K1K2Kps  0 


iqloKps  0 

-Vdr'oKps  0 

Vq'^Kps  0 

0  ~Cg 

KjKjKj,,  0 

0  0 

1  0  . 


(183) 


(185) 


where  the  additional  constants  are  given  by 

dj  =  Cj  +  KpwmKpqC4 

do  =  Kpwm  kjq  ^4 
d3  =  Kpwm  Kpq  C4 
d4  =  Cj  +  KpwmKpdC4 
d5  =  Kpwm  IQd  C4 

aj  =  -c3Vdro  -  i<ko(l-Kps)  -  d3 Kt K2 Kps 
a3  —  C3  Vqro  Iqso  (  1  —  ^ps  ) 
a3  =  -  i&o  Kis  +  d3K!K2Kis 


(186) 

(187) 

(188) 

(189) 

(190) 

(191) 

(192) 

(193) 


The  operating  point  quantities  listed  in  the  above  are  found  once  the  desired  load  torque 
and  commanded  speed  are  specified.  For  instance. 
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(197) 


iqso  —  Vlas.pko  idso 

The  rotor  flux  linkages  per  second  operating  points  are  then  found  from 
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Let’s  consider  a  design  example  for  the  machine  parameters  listed  in  Table  14.  For 
the  given  rated  speed  and  rated  voltage,  the  rated  magnetizing  current  may  be  evaluated  by 
circuit  analysis 

W  =  0.9684  A 

We  also  must  choose  =  1.7  A  and  a  corresponding  limit  on  the  slip  frequency  of 
cosUim  =  30  rad/sec.  These  quantities  are  required  to  calculate  K!  in  the  linearized  system. 
We  also  will  assume  that  Vdc  =  850  V  and  At  =  10  V.  If  we  choose  the  coefficient  of  damping 
to  be  B= 0.001  N-m-sec,  the  steady-state  load  torque  to  be  TLo=0.0N-m,  and  the  steady-state 
desired  speed  to  be  cor*  = 300  rad/sec,  the  following  operating  points  are  identified  by 
calculation 


ie  = 

1qso 

0.2803  A 

idso  = 

0.9684  A 

II 

"g. 

S? 

0.0  V 

II 

> 

144.29  V 

®eo  = 

306.028  rad/sec 

Ias,pko 

=  1.008  A 

Note  that  since  this  is  a  4-pole  machine,  on*  =  300 rad/sec  corresponds  to  a  running 
speed  of  150  rad/sec  or  1432  rpm.  These  values  may  then  be  substituted  along  with  the 
machine  parameters  into  the  linearized  system  and  the  eigenvalues  evaluated  for  various 
gains.  The  following  set  of  gains,  Kpq  =  Kpd  =10,  Kiq  =  Kid  =50000,  Kps  =  5.0,  and  Kis  =  75.0, 
resulted  in  the  following  set  of  eigenvalues: 

-4001  +/-  j  4919 
— 3805  +/-  j4683 
-15.3  +/-  j  5.7 
— 385 
— 16.3 

Note  that  all  of  the  eigenvalues  are  comfortably  in  the  open  left-half  plane  implying  that 
the  selected  operating  point  is  stable  and  should  exhibit  acceptable  transient  response.  The 
first  two  sets  of  eigenvalues  are  largely  influenced  by  the  current  control  gains  while  the 
remaining  four  eigenvalues  are  influenced  by  the  speed  control  parameters  and  operating 
point.  It  should  be  clear  that  the  current  control  dynamics  are  much  faster  than  the 
mechanical  dynamics  and  thus  our  original  assumptions  are  true.  A  decoupled  analysis  of 
the  inner  and  outer  control  loops  may  be  performed  to  arrive  at  a  more  analytical  approach 
to  determining  the  control  gains,  however,  that  will  not  be  presented  in  this  document.  To 
investigate  the  robustness  of  the  design  let’s  consider  another  operating  point,  TLo  =  0.3N-m 
and  co*  =  75.0 rad/sec,  the  new  operating  point  values  are 


0.3504  A 
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i&>  =  0.9684  A 
Vqreo  =  0.0V 
V&  =  144.29  V 
a>eo  =  82.535  rad/sec 

las.pko  =  103  A 

The  corresponding  eigenvalues  are  located  at 

-3852  +/-  j  4837 
-3988  +/-  j4727 
-16.7  +/-  j  8.2 
-311 
-16.0 

and  the  response  is  anticipated  to  be  acceptable. 

The  next  step  is  to  evaluate  the  system  using  a  detailed  ACSL  simulation.  First  a 
couple  subtleties  must  be  pointed  out.  For  the  system  to  develop  the  requisite  torque  to 
begin  accelerating  the  machine,  it  is  advisable  to  first  establish  the  air-gap  flux  by  bringing 
up  ijf  to  its  desired  value.  Next,  since  the  machine  torque  is  limited  by  the  current  bound, 
the  machine  can  only  accelerate  so  quickly.  It  is  generally  wise  to  ramp  up  the  speed  at  a 
rate  that  the  machine  can  logically  follow.  Furthermore,  the  speed  control  integrator 
should  have  anti-windup  incorporated  to  avoid  large  overshoots  and  settling  times. 

An  ACSL  simulation  was  developed  and  the  following  study  conducted.  The  initial 
0.1  seconds  are  used  to  build  up  the  air-gap  flux  in  the  motor  while  it  is  at  rest.  The. 
commanded  speed  is  then  ramped  up  from  Orad/sec  to  300rad/sec  over  the  next  0.6  seconds. 
The  simulation  is  then  run  till  0.9  seconds  to  assess  the  settling  of  the  rotor  speed  and  other 
variables.  In  this  initial  study,  the  load  torque  is  set  to  zero  and  only  friction  is  present. 
The  previous  set  of  parameters  are  assumed  and  a  switching  frequency  of  5  kHz  is 
employed. 

The  following  figures  illustrate  the  dynamics  (Fig.  17-20).  Note  that  the  actual  rotor 
electrical  angular  velocity  tracks  the  desired  value  and  only  experiences  a  slight  overshoot 
(303.7  rad/sec)  at  the  transition  point  (0.7  sec).  It  is  also  clear  that  the  torque  follows  the 
commanded  slip  frequency  as  expected.  In  addition,  the  limit  on  the  slip  frequency 
(30  rad/sec)  translates  to  a  limit  on  the  torque  (1.48  N-m)  and  a  corresponding  limit  on  iqf. 
The  current  control  is  fast  and  accurate  as  is  illustrated  by  the  d-axis  current  being  locked 
around  it’s  desired  value  (0.9684 A).  Harmonics  are  apparent  in  each  of  the  SYN-RF 
currents  and  can  be  further  reduced  by  increasing  the  switching  frequency  from  5  kHz.  The 
next  set  of  curves  illustrate  a  step  change  in  load  torque  from  ON-m  up  to  0.5  N-m.  The 
drive  responds  quickly  and,  following  a  1.5  rad/sec  dip  in  speed,  regains  operation  at  the 
commanded  speed.  Also,  the  air-gap  flux  remains  very  well  regulated  at  the  rated  value. 
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A  second  study  was  then  performed  in  which  the  commanded  speed  was  reduced  to 
75  rad/sec  (Figs.  21-22).  Once  again  the  commanded  speed  was  ramped  up  following  an 
initial  period  in  which  the  air-gap  flux  was  brought  up  to  its  rated  value.  The  machine 
accelerates  under  current-control  up  to  its  new  steady-state  speed  with  little  overshoot.  At 
0.9  seconds,  a  load  torque  of  0.3  N-m  was  stepped  onto  the  machine.  As  illustrated  in  the 
figure,  the  drive  responds  quickly  and  returns  the  motor  to  operation  as  75  rad/sec. 


2.  Indirect  Method  of  Vector  Control 

The  second,  and  final,  implementation  for  an  induction  machine  startup  algorithm 
involves  using  the  indirect  method  of  vector  control.  Vector  control  achieves  near- 
instantaneous  torque  control  by  properly  partitioning  the  stator  current  into  a  flux¬ 
controlling  component  and  a  torque-controlling  component.  The  flux-controlling 
component  is  typically  held  fixed  (unless  field-  weakening  mode  is  required)  so  that  the 
rotor  flux  is  maintained  at  its  optimal  rated  level.  The  torque-controlling  component  of  the 
stator  current  is  adjusted  based  on  the  output  of  the  speed-control  loop.  This  decoupling 
control  is  best  explained  using  the  SYN-RF  representation  of  the  induction  machine.  The 
point  of  the  SYN-RF  transformation,  in  addition  to  eliminating  time-varying  inductances, 
was  to  create  q  and  d-windings  that  were  mutually  orthogonal  so  that  q-currents  only  effect 
q-flux  and  d-currents  only  effect  d-flux.  As  a  consequence,  Our  decoupling  of  torque  and 
flux  control  is  realized  by  recognizing  that  we  need  either  the  q  or  d-component  of  the  flux 
to  be  zero  so  that  the  corresponding  current  component  will  not  effect  the  flux.  It  has  been 
shown  that  instantaneous  torque  control  is  achieved  by  implementing  the  control  with 
regards  to  the  rotor  flux  instead  of  either  the  stator  flux  or  the  air-gap  flux.  Thus,  we 
impose 

\|/qre  =  0  (200) 

and 


(201) 


Therefore  iqes  is  used  to  control  the  torque  and  i<£  is  used  to  set  the  rotor  flux.  There  are 
two  approaches  to  realizing  these  constraints:  the  direct  and  the  indirect  method.  In  the 
direct  method,  the  spatial  position  of  the  peak  rotor  flux  is  located.  This  is  where  we  must 
place  the  SYN-RF  d-axis.  Thus,  we  arrive  at  a  0e*  which  enables  us  to  form  iqes  and  ides  from 
the  measured  stator  currents  and  then  also  transform  the  SYN-RF  control  signals  back  to 
abc-quantities  for  performing  STPWM.  In  the  indirect  method,  a  necessary  and  sufficient 
condition  between  iqes,  i<£,  and  <»e-a>r  is  established  from  the  modeling  equations.  This 
results  in  a  desired  slip  angle  to  which  we  add  the  actual  rotor  electrical  position  to  derive 
the  required  SYN-RF  angle.  The  0*  is  then  used  as  described  previously. 
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The  slip  relationship  that  must  hold  for  the  q-axis  rotor  flux  to  be  held  at  zero  is 

rrXcOb  iqf  (202) 

coe  -  cor  =  — - 7  v  ' 

X„  Vdr 

If  we  choose  to  hold  idf  constant  so  that  =  Xmidf,  then  it  follows  that  the  required 
SYN-RF  angle  is  given  by 


=  0r  +  J 


rrC0b 


_3L 

:  e* 
Ids 


dC 


(203) 


Therefore,  the  principle  difference  between  the  vector  control  implementation  and  the 
previous  approach  is  now  the  speed  control  loop  output  is  iqf  and  co*  is  no  longer  required 
since  we  are  directly  forming  6e*.  The  nonlinear  mapping  between  (os,  and  1^  is  no  longer 
needed  as  is  the  partioning  of  1^  into  iqf  and  idf . 

The  attractice  features  of  vector  control  are  realized  as  long  as  the  control  remains 
tuned,  that  is,  as  long  as  the  parameters  used  in  the  slip  relationship  closely  match  those  of 
the  machine.  Typically,  adaptive  algorithms,  such  as  an  extended  Kalman  filter,  may  be 
merged  with  the  vector  control  to  provide  updates  on  the  rotor  resistance  which  is  a  very 
temperature  dependent  parameter.  Such  an  expense  is  not  warranted  for  drives  which  do 
not  require  instantaneous  torque  response  and  is  not  considered  here. 

The  design  of  the  current  control  gains  and  speed  control  gains  may  be  accomplished 
in  a  manner  analogous  to  the  previous  section.  First,  note  that  the  equations  for  the 
induction  machine  will  be  the  same.  Next,  the  slip  frequency  constraint  is  forcing 


o>e 


rr  «>b 

x^r 


*qs 


+  C0r 


=  T,  iqes 


+  C0r 


(204) 


and  thus  linearizing 


Acoe*  =  T,  A  iqf  +  A(0r 


(205) 


The  output  of  the  speed  control  loop  is  similar  to  before 

iqf  =  Kps  (to*  -  ©r)  +  KisXs 

or  linearizing 


(206) 


A  iqf  =  KpS(Acor*-A(Or)  +  KisAxs 


(207) 


Inserting  the  current  control  dynamics  and  the  linearization  of  the  induction  machine 
dynamics  results  in 


-^j-A£im2  =  Aim2  A5tim2  +  Bim2  Atfim2 


(208) 


Aim2  — 


(209) 
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Mm2  =  [Ac°r* 


(212) 


where  the  additional  constants  are 

S4  =  — idso (  1  —  TjKps)  —  C3  Vdro  —  ^3  Kps 
=  lqso(l  OTj  ICpS )  +  C3  ^|!qr0 
^6  =  -i<£oTiKjS  +  d3Kis 


(213) 

(214)  • 

(215) 


In  order  to  determine  the  operating  points  to  be  inserted  in  the  above  linearized 
system,  the  peak  rotor-referred  flux  must  be  identified  by  analyzing  the  steady-state 
equivalent  circuit  for  rated  conditions.  This  yields 

xfao  =  143.725  V 

Also  if  the  control  is  tuned,  the  necessary  and  sufficient  condition  guarantees  that 


Vqn 


=  ov 


The  d-axis  stator  current  is  fixed  by  the  desired  rotor  flux  according  to 

Vdro 


idso  ~ 


X„ 


=  0.9646  A 


The  q-axis  stator  current  is  established  from  the  cummulative  load  on  the  shaft, 
during  tuned  vector  control 


Since 


Te  = 


3  P  Xm  iqSO 
4©bXr;  Vdro 


—  Tlo  +  B  B)r0 


(216) 


the  q-axis  stator  current  is  found  readily.  Let’s  consider  a  numerical  example  and  illustrate 


-62- 


some  effective  gains.  Assume  that  =  ON-m  and  <o*  =  300rad/sec.  This  leads  to 

iqeso  =  0.2814  A 
(Oeo  =  306.075  rad/sec 

Substituting  into  the  linearized  system  equations,  gains  may  now  be  selected  and  the 
eigenvalues  evaluated.  Once  again,  it  was  assumed  that  vdc  =  850V  and  At=10V.  The 
following  set  of  gains,  Kpq  =  Kpd=lO,  Kiq  =  Kid  =50000,  Kps  =  0.1,  and  Kis  =  5.0,  resulted  in  the 

following  set  of  eigenvalues: 

-4102  +/-  j4890 
-3853  +/-  j  4587 
— 44.7  +/-  j49.5 
-20.8  +/-  j  6.1 

A  variety  of  operating  points  were  then  assessed  and  it  was  found  that  the  eigenvalues 
were  largely  fixed  in  the  complex  plane,  with  the  dominant  eigenvalues  at  around  -20. 

The  vector  control  was  next  to  be  assessed  via  detailed  ACSL  simulation.  Similar 
provisions  were  placed  on  the  speed  control  to  ensure  against  integrator  anti-windup. 
Initially,  0.1  seconds  are  reserved  to  allow  the  rotor  flux  to  be  established  at  which  point  the 
commanded  speed  is  ramped  up  over  a  0.6  second  interval  and  then  held  fixed  for  an 
additional  0.2  seconds.  The  results  are  documented  in  the  following  figures. 

In  the  study  (Figs.  23-26),  a  0.1  second  time  period  is  used  to  build  up  the  rotor  flux  to 
its  prescribed  value.  The  commanded  speed  is  then  ramped  up  to  300  rad/sec  over  a  0.6  sec 
interval.  The  study  is  run  till  t= 0.9  sec.  The  drive  responds  with  excellent  performance 
only  a  slight  2%  overshoot.  At  0.9 sec  the  load  torque  is  stepped  from  ON-m  up  to  0.5 N-m. 
The  eigenvalue  selections  result  in  a  slightly  overdamped  response  which  settles  out  in 
0.2  sec.  The  transient  in  the  rotor  speed  is  quite  acceptable.  The  change  is  rotor  flux  is 
hardly  noticeable. 
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Figure  19.  Transient  Behavior  During  a  Load  Torque  Change 


Figure  20.  Current  and  Flux  Waveforms  Corresponding  to  Figure  19. 
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Figure  25.  Transient  Behavior  During  a  Load  Torque  Change 


~0. 90000 


0.93333 


0.96667  1.00000 

T 


1.03333 


1.06667 


1. 10000 


Figure  26 •  Current  and  Flux  Waveforms  Corresponding  to  Figure  25. 
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C.  Concluding  Comments 

The  Navy  is  actively  developing  Power  Electronic  Building  Blocks  (PEBBs)  in 
support  of  its  initiatives  into  DC  Zonal  Electric  Distribution  (DC  ZEDS).  DC  ZEDS,  in 
turn,  is  a  component  of  the  larger  Integrated  Power  System  (IPS)  initiative  that  seeks  to 
power  both  ship  service  and  ship  propulsion  from  a  common  set  of  prime  movers.  The 
principle  PEBB-like  components  under  development  currently  are  the  Ship  Service 
Converter  Module  (SSCM)  and  the  Ship  Service  Inverter  Module  (SSIM)  as  specified  in 
DC  ZEDS.  These  devices  are  basically  fast-switching,  high-bandwidth,  high-power  buck 
choppers  and  Pulse-Width-Modulation  (PWM)  inverters,  respectively.  Currently,  the 
Auxiliary  Resonant  Commutated  Pole  (ARCP)  inverter  is  the  topology  of  choice  due  to  its 
inherent  soft-switching  advantages  while  a  Zero- Voltage-Switching  (ZVS)  buck  chopper  is 
being  investigated  to  replace  its  hard-switched  counterpart. 

The  purpose  of  this  report  was  to  document  several  algorithms  which  may  be 
employed  to  control  a  conventional  three-phase  PWM  inverter.  It  is  assumed  that  the 
soft-switching  aspect  of  the  ARCP  inverter  is  being  handled  internal  to  the  unit.  The 
report  first  documented  two  approaches  to  developing  the  actual  modulation  signals.  First, 
the  discretization  of  Sine-Triangle  PWM  (STPWM)  was  considered  and  the  issues 
regarding  implementation  with  the  PEBB  Universal  Controller  Generation  0  (PUC-GENO) 
were  highlighted.  Specific  coding  details  are  documented  in  both  Hanson  and  Floodeen. 
Next,  Space  Vector  Modulation  (SVM)  was  outlined  and  various  equations  required  for 
implementation  identified.  SVM  has  an  advantage  of  being  able  to  realize  slightly  larger 
voltages  in  the  linear  modulation  range  but  also  requires  far  more  calculation  capacity  or 
more  extensive  look-up  tables.  This  poses  a  trade-off.  For  very  complex  feedback  control 
algorithms,  the  extra  overhead  involved  in  SVM  may  preclude  its  selection. 

The  next  section  of  the  report  documented  current  and  voltage  control  strategies  for  a 
three-phase  PWM  inverter.  In  particular  with  regards  to  the  current  control,  algorithms 
were  presented  which  were  implemented  in  both  the  Stationary  Reference  Frame  (STA-RF) 
and  the  Synchronous  Reference  Frame  (SYN-RF).  For  generality,  a  standard  three-phase 
resistive-inductive  (RL)  load  was  assumed  for  convenience.  The  STA-FR  control  is  more 
easily  implemented  since  there  are  no  angle-dependent  transformations  required;  however, 
steady-state  accuracy  is  very  parameter  dependent  and  therefore  care  must  be  exercised. 
Modeling  equations  were  derived  for  the  STA-RF  control  and  design  example  used  to 
illustrate  the  selection  of  the  feedback  gains.  A  detailed  ACSL  simulation  was  then  used 
to  illustrate  the  system  dynamics.  The  SYN-RF  control  is  slightly  more  complicated  due 
to  the  fact  that  the  feedback  signals  must  be  transformed  to  the  SYN-RF  via  a 
diffeomorphic  angle-dependent  transformation.  In  addition,  signals  internal  to  the  control 
mus  then  be  mapped  back  to  abc-quantities  for  STPWM  or  to  STA-RF-quantities  for  SVM. 
However,  a  SYN-RF  control  was  successfully  implemented  in  the  NPS  Power  Systems 
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Laboratory  using  the  dSPACE  control  development  system  (Frasz,  1998)  and  its  advantage 
far  outweighs  the  added  complexities.  In  particular,  the  integral-action  in  the  SYN-RF 
guarantees  zero  steady-state  error  since  sinusoidal  quantities  appear  as  constants  in  the 
SYN-RF.  Modeling  equations  were  derived  for  the  SYN-RF  control  and  a  similar  design 
example  used  to  illustrate  the  selection  of  feedback  gains.  An  ACSL  simulation  was 
implemented  and  used  to  study  the  dynamic  performance.  Results  were  found  to  be 
acceptable. 

Most  shipboard  applications  will  probably  require  a  well-regulated  output  voltage. 
The  next  section  of  the  report  documented  a  sophisticated  inverter  output  voltage  control 
which  relies  heavily  on  reference  frame  theory.  It  incorporates  an  outer  voltage  control 
loop  and  an  inner  current  control  loop  and  various  feedforward  signals.  The  modeling 
equations  were  set  forth,  a  design  example  used  to  illustrate  the  selection  of  feedback 
gains,  and  a  detailed  ACSL  simulation  used  to  document  the  dynamic  performance.  This 
algorithms  requires  a  number  of  feedback  signals  but  results  in  superior  dynamic 
performance  as  illustrated  in  the  aforementioned  studies.  A  systematic  algorithm  for 
determining  the  feedback  gains  was  set  forth  and  shown  to  be  independent  of  the  load 
parameters  themselves. 

The  final  section  of  the  report  dealt  with  two  algorithms  for  handling  the  startup  of  a 
three-phase  induction  machine  load.  Both  algorithms  utilized  an  outer  speed  control  loop 
and  an  inner  current  control  loop.  The  first  algorithm  used  the  output  of  the  outer  speed 
control  loop  to  program  a  desired  slip  frequency  which,  in  turn,  was  bounded  to  limit  the 
torque  and  current.  The  commanded  slip  frequency  was  added  to  the  rotor  electrical 
angular  velocity  to  arrive  at  the  desired  inverter  switching  frequency  (in  rad/sec).  In 
addition,  the  commanded  slip  frequency  signal  was  then  used  to  map  into  a  commanded 
stator  current  amplitude  in  order  to  achieve  constant  air-gap  flux  in  the  steady  state.  The 
current  amplitude  was  then  apportioned  into  commanded  synchronous  reference  frame 
currents  and  sent  to  the  inner  current  control,  which  was  simply  the  SYN-RF  current 
control  from  a  previous  section.  The  modeling  equations  were  set  forth  and  a  linearized 
system  derived  to  aid  in  the  selection  of  control  gains  for  a  representative  machine  from 
the  NPS  Power  Systems  Laboratory.  A  detailed  ACSL  simulation  was  formulated  to 
illustrate  the  dynamics. 

Finally,  an  alternative  startup  strategy  was  derived  using  the  indirect  method  of  vector 
control.  The  theory  was  introduced,  the  equations  assembled,  a  linearized  system 
formulated,  a  design  example  proffered,  and  an  ACSL  simulation  generated  and  exercised. 
The  basic  difference  between  this  algorithm  and  the  previous  is  that  a  necessary-and- 
sufficient  condition  on  the  slip  frequency  is  derived  which  guarantees  near-instantaneous 
torque  response.  This  condition  is  used  in  a  feedforward  sense  together  with  a  measure  of 
the  rotor  position  to  derive  a  desired  SYN-RF  angle.  This  angle  is  then  used  directly  in 
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the  control  to  map  the  SYN-RF  signals  back  to  abc-modulating  signals.  The  simulation 
study  results  illustrated  acceptable  performance. 

PWM  Inverter  control  algorithms  must  be  derived  with  vision  of  all  anticipated 
loading  conditions.  With  this  in  mind,  additional  considerations  must  be  kept  in  mind.  In 
particular,  paralleling  of  inverter  modules  to  share  a  load  proportionately  is  a  key  issue. 
Also,  the  impact  of  input  filters  on  overall  system  stability  must  be  an  issue  which  is 
heeded. 
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